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I .  OBJECTIVES 

The  objectives  of  this  project  were  to  investigate  the  full-wave 
(Maxwellian)  analysis  for  solving  and  computing  the  electromagnetic 
behavior  of  microwave  integrated  circuits. 

II.  SUMMARY  OF  WORK  DONE 

The  majority  of  the  work  was  done  in  the  first  two  years  of  this 
project.  Summary  reports  for  the  first  year's  work  and  the  second  year's 
work  are  included  to  give  details.  The  final  six  month's  period  was  used 
primarily  to  prepare  final  manuscripts  for  submission  to  journals,  and  to 
give  a  preliminary  organization  for  Mr.  Hsu's  Ph.D.  dissertation.  All 
publications  to  date  and  the  preliminary  version  of  Mr.  Hsu's  dissertation 
are  appended  to  this  report.  The  following  is  a  list  of  these  appendices: 

1.  Annual  Report  for  the  period  January  1,  1988  to  December  31,  1988. 

2.  Annual  Report  for  the  period  January  1,  1989  to  December  31,  1989. 

3.  "Quasi-Static  Analysis  of  a  Micrdstrip  Via  Through  a  Hole  in  a 
Ground  Plane,"  by  Taoyun  Wang,  R.  F.  Harrington,  and  J.  R.  Mautz, 

IEEE  Trans,  on  MTT,  Vol.  36,  No.  6,  pp.  1008-1013,  June  1988. 

4.  "The  Inductance  Matrix  of  a  Multiconductor  Transmission  Line  in 
Multiple  Magnetic  Media,"  J.  R.  Mautz,  R.  F.  Harrington,  and  C-I  G.  Hsu, 

IEEE  Trans,  on  MTT,  Vol.  36,  No.  8,  pp.  1293-1295,  August  1988. 

5.  "Boundary  Integral  Formulations  for  Homogeneous  Material  Bodies," 

R.  F.  Harrington,  Journ.  of  EM  Waves  and  Appl,  Vol.  3,  No.  1,  pp.  1-5,  1989. 

6.  The  Excess  Capacitance  of  a  Microstrip  Via  in  a  Dielectric 
Substrate,"  IEEE  Trans,  on  CAD,  Vol.  9,  No.  1,  pp.  48-56,  January  1990. 

7.  "A  Stable  Integral  Equation  for  Electromagnetic  Scattering  from 
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Homogeneous  Dielectric  Bodies,"  J.  R.  Mautz,  IEEE  Trans,  on  AP,  Vol.  37, 
No.  8,  pp.  1070-1071,  August  1989. 

8 .  "On  the  Location  of  Leaky  Nave  Poles  for  a  Grounded  Dielectric 
Slab,"  C-I  G.  Hsu,  R.  F.  Harrington,  J.  R.  Mautz,  and  T.  K.  Sarkar, 
submitted  to  IEEE  Trans,  on  MTT . 

9.  "Fundamental  Mode  Dispersion  of  a  Microstrip  Transmission  Line," 

G-I  G.  Hsu,  R.  F.  Harrington,  J.  R.  Maut  -C  Jan,  submitted  to 

IEEE  Trans,  on  MTT. 

10.  "On  the  Dynamic  Model  for  Passive  Microwave  Printed  Circuits," 
C-I  G.  Hsu,  dissertation  proposal,  June  1990. 

III.  PERSONNEL 

The  following  persons  have  worked  on  this  project : 

Roger  F.  Harrington,  Principal  Investigator,  Professor  of  Electrical 
Engineering,  Syracuse  University. 

Joseph  R.  Mautz,  Research  Engineer,  Department  of  Electrical  and 
Computer  Engineering,  Syracuse  University. 

Chung-I  Gavin  Hsu,  Research  Assistant,  Graduate  Student,  Department 
of  Electrical  and  Computer  Engineering,  Syracuse  University. 

Ing-Chieh  Jan,  Research  Assistant,  Graduate  Student,  Department 
of  Electrical  and  Computer  Engineering,  Syracuse  University. 
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I.  OBJECTIVES 


The  objectives  of  this  project  are  to  obtain  full-wave  (Maxwellian) 
analyses  suitable  for  computing  the  electromagnetic  behavior  of  micro- 
wave  integrated  circuits.  Initial  work  has  been  on  (a)  attempting  to  ex¬ 
tend  the  quasi-static  analyses  to  include  correction  terms  to  allow 
application  at  higher  frequencies,  and  (b)  full-wave  analyses  to  two- 
dimensional  problems  such  as  multiconductor  transmission  lines  in  multi¬ 
ple  dielectric  media. 

II.  SUMMARY  OF  FIRST  YEAR*  S  WORK 


A)  Preliminary  to  extending  the  analyses  to  higher  frequencies,  we 
completed  some  quasi-static  solutions.  The  following  topics  have  led  to 
papers  published  or  pending,  as  listed  in  Section  IV. 

(1)  Quasi-static  analysis  of  a  microstrip  via  through  a  hole 
in  a  ground  plane  and  in  a  dielectric  substrate.  Abstract  of  [1]:  The 
equivalent  circuit  of  a  via  which  connects  two  semi-inf initely  long  trans¬ 
mission  lines  through  a  circular  hole  in  a  ground  plane  is  considered. 

The  iT-type  equivalent  circuit  consists  of  two  excess  capacitances  and  an 
excess  inductance.  They  are  quasi-static  quantities  and  thus  are  computed 
statically  by  the  method  of  moments  from  the  integral  equations.  The 
integral  equations  are  established  by  introducing  a  sheet  of  magnetic  cur¬ 
rent  in  the  electrostatic  case  and  a  layer  of  magnetic  charge  in  the  mag¬ 
netostatic  case.  Parametric  plots  of  the  excess  capacitances,  the  excess 
inductance,  and  the  characteristic  admittance  of  the  via  are  given  for 
reference. 

Abstract  of  [4]:  The  equivalent  circuit  of  a  via  which  connects  two 
semi- inf initely  long  microstrip  transmission  lines  imbedded  in  a  dielec¬ 
tric  medium  above  a  ground  plane  is  considered.  The  T-type  equivalent 
circuit  consists  of  an  excess  capacitance  and  an  excess  inductance.  The 
excess  inductance  of  the  via  is  the  same  as  the  one  computed  in  the  case 
of  free  space  (without  a  substrate).  The  excess  capacitance  of  the  via  is 
computed  quasi-statically  by  the  method  of  moments  and  the  image  method 
from  the  integral  equations.  It  converges  rapidly  as  the  number  of  image 
terms  is  increased.  Parametric  plots  of  the  excess  capacitance  of  the  via 
are  given  for  reference. 

(2)  The  equivalent  circuit  of  a  microstrip  crossover  with  and 
without  a  dielectric  substrate.  Abstract  of  [4]:  The  equivalent  circuit 
of  a  raicrostrip  crossover  is  found.  Integral  equations  are  obtained  for 
the  densities  of  excess  charge,  and  these  equations  are  solved  by  the  method 
of  moments.  Introduction  of  a  specified  transverse  distribution  of  charge, 
which  satisfies  the  edge  condition,  reduces  the  computing  time  dramatically 
while  the  accuracy  remains  excellent.  Several  plots  of  the  excess  charge 
densities  are  provided  along  with  numerical  values  of  lumped  excess  capaci¬ 
tances  . 

Abstract  of  [5]:  A  quasi-static  analysis  is  carried  out  to  examine 
the  capacitive  coupling  between  two  nonintersecting  microstrip  lines  above 
a  ground  plane  and  in  a  dielectric  substrate.  The  charge  density  along 


5 


the  width  of  each  strip  is  described  using  a  prescribed  charge  distribu¬ 
tion.  A  pair  of  coupled  integral  equations  is  derived  and  solved  via 
the  method  of  moments  to  obtain  the  excess  charge  densities  .  The  lumped 
excess  capacitances  are  computed  and  compared  to  the  ones  obtained  using 
wire  lines  with  radii  equal  to  the  equivalent  radii  of  the  strips. 

(3)  The  inductance  matrix  of  a  multiconductor  transmission 
line  in  multiple  magnetic  media.  Abstract  of  [2]:  Consider  a  multicon¬ 
ductor  transmission  line  consisting  of  Nc  conducting  cylinders  in  inhomo¬ 
geneous  media  consisting  of  Nj  homogeneous  regions  with  permeabilities  Uj_ 
and  permittivities  The  inductance  matrix  [L]  for  the  line  is  obtained 

by  solving  the  magnetostatic  problem  of  Nc  conductors  in  Nd  regions  with 
permeabilities  u^.  The  capacitance  matrix  [ C ]  for  the  line  is  obtained 
by  solving  the  electrostatic  problem  of  Nc  conductors  in  Nd  regions  with 
permittivities  £^.  It  is  shown  that  [L]  *  UQ£0  where  [ C ' ]  is  the 

capacitance  matrix  of  an  auxiliary  electrostatic  problem  of  Nc  conductors 
in  Nd  regions  with  relative  permittivities  set  equal  to  the  reciprocals 
of  the  relative  permeabilities  of  the  magnetostatic  problem,  i.e.,  e!/e  = 

U  /u .  • 
o  i 

B)  Work  has  been  done  on  the  following  full-wave  solutions: 

(1)  The  full-wave  solution  for  the  microstrip  transmission  line 
has  been  investigated  using  the  volume  integral  formulation.  This  gives  an 
integral  equation  formulation  where  the  unknowns  are  the  conduction  current 
and  charge  on  conducting  boundaries,  the  bound  charge  of  the  dielectric 
boundaries,  and  the  polarization  current  in  the  dielectric  volume.  For 
waves,  the  propagation  constant  is  obtained  from  the  solution  of  a  nonlinear 
eigenvalue  equation.  A  numerical  solution  to  this  equation  has  been  obtained 
to  allow  computation  of  the  propagation  constant,  or,  equivalently,  the 
effective  dielectric  constant.  The  solution  is  complicated,  and  the  computa¬ 
tions  are  time  consuming.  We  are  presently  looking  for  ways  to  simplify  this 
solution. 


(2)  The  full-wave  solution  for  the  microstrip  transmission  line 
has  been  investigated  using  a  surface  integral  formulation.  This  gives  an 
integral  equation  where  the  unknowns  are  the  conduction  current  on  the  con¬ 
ducting  boundaries  and  the  equivalent  electric  and  magnetic  surface  cur¬ 
rents  on  the  dielectric  boundaries.  Again  the  propagation  constant  is 
obtained  from  the  solution  to  a  nonlinear  eigenvalue  equation.  The  numeri¬ 
cal  solution  to  this  equation  has  not  yet  been  successfully  obtained.  Work 
is  continuing  on  this  solution. 

There  are  many  different  surface  integral  formulations  for  any  par¬ 
ticular  problem.  These  are  summarized  in  the  paper  on  boundary  integral 
formulations  for  homogeneous  material  t  lies  ([3]  of  Section  IV).  The 
following  is  an  abstract  of  [3]:  There  are  many  boundary  integral  formu¬ 
lations  for  the  problem  of  electromagnetic  scattering  from  the  transmis¬ 
sion  into  a  homogeneous  material  body.  The  only  formulations  which  give  a 
unique  solution  at  all  frequencies  are  those  which  involve  both  electric 
and  magnetic  equivalent  currents;  and  satisfy  boundary  conditions  on  both 
tangential  E  and  tangential  8.  Formulations  which  involve  only  electric 
(or  magnetic)  equivalent  currents,  and  those  which  involve  boundary  condi- 
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tions  on  only  tangential  E  (or  tangential  H)  are  singular  at  frequencies 
corresponding  to  the  resonant  frequencies  of  a  resonator  formed  by  a 
perfect  conductor  covering  the  surface  of  the  body  and  filled  with  the 
material  exterior  to  the  body  in  the  original  problem. 

I IT.  PLANS  FOR  THE  SECOND  YEAR'S  WORK 

A)  Study  of  the  volume  integral  formulation  for  the  full-wave  solution 
of  2-D  problems  will  continue.  Methods  of  reducing  the  complexity  of  the 
solution  will  be  considered.  In  particular,  a  perturbation  approach,  start¬ 
ing  from  the  quasi-TEM  solution  and  adding  correction  terms  will  be  studied 
and  compared  to  the  direct  solution. 

B)  Study  of  the  surface  integral  formulation  for  the  full-wave  solu¬ 
tion  of  2-D  problems  will  continue.  Again  methods  of  redurng  the  complexity 
of  the  solution  will  be  considered.  However,  this  formula  n  does  not 
appear  to  be  as  well  suited  to  the  perturbation  approach  as  does  the  volume 
integral  formulation  of  A  above. 

C)  Study  of  the  integral  formulations  for  the  full-wave  solution  of 
3-D  problems  will  be  started.  The  particular  approach  to  be  taken  will 
depend  on  the  results  of  A  and  B  above.  It  is  hoped  that  a  3-D  solution 
will  lead  to  a  formulation  useful  for  actual  printed  circuits. 

IV.  PUBLICATIONS 


[1]  T.  Wang,  R.  F.  Harrington,  and  J.  R.  Mautz,  "Quasi-Static  Analysis  of 
a  Microstrip  Via  through  a  Hole  in  a  Ground  Plane,"  IEEE  Trans. , 

vol .  MTT-36 ,  No.  6,  June  1988. 

[2]  J.  R.  Mautz,  R.  F.  Harrington,  and  C-I  G.  Hsu,  "The  Inductance  Matrix 
of  a  Multiconductor  Transmission  Line  in  Multiple  Magnetic  Media," 
IEEE  Trans.,  vol.  MTT-36,  No.  8,  pp.  1293-1295,  Aug.  1988. 

Accepted  for  publication: 

[3]  R.  F.  Harrington,  "Boundary  Integral  Equations  for  Homogeneous 
Material  Bodies,"  Journ.  Electromagnetic  Waves  and  Applications, 
in  press,  to  appear  1988. 

Submitted  for  publication: 

[4]  T.  Wang,  J.  R.  Mautz,  and  R.  F.  Harrington,  "The  Excess  Capacitance 
of  a  Microstrip  Via  in  a  Dielectric  Substrate,"  submitted  to  IEEE 
Transactions  on  Integrated  Circuits  and  Systems. 

[5]  S.  Papatheodorou,  R.  F.  Harrington,  and  J.  R.  Mautz,  "The  Equivalent 
Circuit  of  a  Microstrip  Crossover,"  submitted  to  Transactions  of  the 
Society  for  Computer  Simulation. 


(6]  S.  Papatheodorou,  R.  F.  Harrington,  and  J.  R.  Mautz,  "The  Equivalent 
Circuit  of  a  Microstrip  Crossover  in  a  Dielectric  Substrate,"  sub¬ 
mitted  to  IEEE  Transactions  on  Microwave  Theory  and  Techniques. 


Report : 

[7]  S.  Papatheodorou,  J.  R.  Mautz,  and  R.  F.  Harrington,  "The  Equivalent 
Circuit  of  a  Microstrip  Crossover  in  a  Dielectric  Substrate,"  Report 
TR-88-11,  Department  of  Electdical  and  Computer  Engineering,  Syracuse 
University,  August  1988. 

V.  PERSONNEL 


The  following  is  the  project  staff  for  the  first  year: 

Principal  Investigator:  Dr.  Roger  F.  Harrington, 

Professor  of  Electrical  Engineering. 

Research  Engineer:  Dr.  Joseph  R.  Mautz, 

Department  of  Electrical  ind  Computer  Engineering. 

Research  Assistant:  Mr.  Chung-I  Gavin  Hsu,  Graduate  Student, 

Department  of  Electrical  and  Computer  Engineering. 
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1.  OBJECTIVES 


The  objectives  of  this  project  are  to  obtain  full-wave  (Maxwellian) 
analyses  suitable  for  computing  the  electromagnetic  behavior  of  microwave 
integrated  circuits.  Initial  work  has  been  on  (a)  attempting  to  extend 
the  quasi-static  analyses  to  include  correction  terms  to  allow 
application  at  higher  frequencies,  and  (b)  full-wave  analyses  to  two- 
dimensional  problems  such  as  multiconductor  transmission  lines  in  multi¬ 
ple  dielectric  media.  More  recent  work  has  been  on  (c)  the  solution  of 
three-dimensional  problems . 

II.  SUMMARY  OF  THE  SECOND  YEAR'S  WORK 

A)  Publication  of  the  first  year's  work  as  journal  articles:  The 
following  papers  have  been  published  during  the  second  year: 

1.  Roger  F .  Harrington,  "Boundary  Integral  Formulations  for 
Homogeneous  Material  Bodies,"  Journal  of  Electromagnetic  Waves  and 
Applications,  vol.  3,  No.  1,  pp.  1-15,  1989. 

Abstract  -  There  are  many  boundary  integral  formulations  for  the 
problem  of  electromagnetic  scattering  from  and  transmission  into  a 
homogeneous  material  body.  The  only  formulations  which  give  a  unique 
solution  at  all  frequencies  are  those  which  involve  both  electric  and 
magnetic  equivalent  currents,  and  satisfy  boundary  conditions  on  both  tan¬ 
gential  E  and  tangential  H.  Formulations  which  involve  only  electric  (or 
magnetic)  equivalent  currents,  and  those  which  involve  boundary 
conditions  on  only  tantential  E  (or  tangential  H)  are  singular  at 
frequencies  corresponding  to  the  resonant,  frequencies  of  a  resonator 
formed  by  a  perfect  conductor  covering  the  surface  of  the  body  and  filled 
with  the  material  exterior  to  the  body  in  the  original  problem. 


2.  Joseph  R.  Mautz,  "A  Stable  Integral  Equation  for  Scattering  from 
Homogeneous  Dielectric  Bodies/"  IEEE  Transactions  on  Antennas  and 
Propagation,  vol.  37,  No.  8,  pp.  1070-1071,  August  1989. 

Abstract  -  It  is  shown  that  a  previously  derived  integral  equation 
for  electromagnetic  scattering  from  a  homogeneous  dielectric  body  does 
not  have  a  unique  solution  at  resonant  frequencies  of  the  cavity  formed 
by  making  the  surface  S  of  th-  body  perfectly  conducting  and  filling  the 
region  internal  to  S  with  the  external  medium.  This  integral  equation 
was  formulated  so  that  an  equivalent  electric  current  radiates  in  the 
presence  of  the  homogeneous  external  medium  to  produce  the  scattered 
field  external  to  the  body.  A  combination  of  equivalent  electric  and 
magnetic  currents  is  used  to  formulate  an  integral  equation  whose 
solution  is  always  unique. 

3.  Taoyun  Wang,  J.  R.  Mautz,  and  R.  F.  Harrington,  "The  Excess 
» 

Capacitance  of  a  Microstrip  Via  in  a  Dielectric  Substrate,"  IEEE 
Transactions  on  Computer-Aided  Design,  vol.  9,  No.  1,  pp.  48-56,  January 
1989. 

Abstract  -  The  equivalent  circuit  of  a  via  which  connects  two  semi- 
infinitely  long  microstrip  transmission  lines  imbedded  in  a  dielectric 
medium  above  a  ground  plane  i3  considered.  The  T-type  equivalent  circuit 
consists  of  an  excess  capacitance  and  an  excess  inductance.  The  excess 
inductance  of  the  via  is  the  same  as  the  one  computed  in  the  case  of  free 
space  (without  a  substrate! .  The  excess  capacitance  of  the  via  is 
computed  quasi-statically  by  the  method  of  moments  and  the  image  method 
from  the  integral  equations.  It  converges  rapidly  as  the  number  of  image 
terms  is  increased.  Parametric  plots  of  the  excess  capacitance  of  the 
via  are  given  for  reference . 


B)  Two-dimensional  full-wave  solutions  of  printed  circuits:  The 
surface  integral  formulation  was  found  to  be  more  suitable  for 
computation  than  was  the  volume  integral  formulation.  Hence,  this 
solution  was  used  for  the  final  computations.  A  Technical  Report  was 
written:  "Fundamental  Mode  Dispersion  of  a  Microstrip  Transmission 
Line,"  by  Chung-I  G.  Hsu,  R.  F.  Harrington,  J.  R.  Mautz,  and  I-C  Jan, 
Report  TR-89-6,  November  1989.  A  paper  by  the  same  name  and  authors  has 
been  submitted  to  the  IEEE  Transactions  on  Microwave  Theory  and 
Techniques . 

Abstract  -  A  spectral  domain  type  solution  with  nonumiform  subdomain 
basis  functions  is  implemented  to  investigate  the  dispersion  character¬ 
istics  of  the  fundamental  mode  of  a  microstrip  transmission  line.  Agree¬ 
ment  with  the  data  available  in  the  literature  is  demonstrated  in  several 

plots,  where  only  the  longitudinal  current  is  used  in  numerical 
% 

computation.  This  numerical  scheme  proves  to  be  an  alternative,  or 
rather  a  supplement,  to  various  well-developed  methods  in  this  area. 

C)  Integral  formulations  and  solutions  for  the  three-dimensional 
printed  circuit  problem.  We  have  obtained  a  formulation  for  the  3-D 
circuit  as  an  extension  of  the  arbitrarily-shaped  wire  formulation  for  the 
homogeneous  dielectric  case.  (For  example,  see  Chapter  4  of  FIELD 
COMPUTATION  BY  MOMENT  METHODS,  R.  F.  Harrington,  Krieger  Publishing  Co., 
1982)  .  The  solution  for  the  printed  circuit  problem  is  much  more  diffi¬ 
cult  than  that  for  the  homogeneous  dielectric  problem,  because  now  the 
Green's  function  contains  a  Sommerfeld  integral.  Some  approximate  and 
numerical  methods  for  evaluating  the  integral  have  been  obtained  in  the 
literature,  but  the  solution  is  still  complicated.  We  intend  to  use  an 


"equivalent  radius"  simplification  to  obtain  a  wire-type  solution  for 


arbitrary  printed  circuits.  The  equivalent  radius  will  be  obtained  from 
the  simpler  2-D  problem  of  a  printed  circuit  transmission  line. 

D)  Software  packages  for  printed  circuit  multiconductor  transmission 
lines:  Three  menu-driven,  user-friendly  software  packages  have  been  pub¬ 
lished  by  Artech  House.  These  are: 

(1)  "Matrix  Parameters  of  Multiconductor  Transmission  Lines,"  by  A. 
A.  R.  D jordjevic,  R.  F.  Harrington,  T.  K.  Sarkar,  and  M.  B.  Bazdar.  This 
program  calculates  the  quasistatic  matrix  parameters  of  multiconductor 
transmission  lines,  i.e.,  matrices  of  inductance,  capacitance,  resist¬ 
ance,  and  conductance.  The  program  can  treat  multiple  dielectric  layers 
as  well  as  multiple  conductors  of  arbitrary  shape.  A  complete  user's 
manual  accompanies  the  computer  program. 

(2)  "Time-Domain  Response  of  Multiconductor  Transmission  Lines,"  by 

A.  R.  D jordjevic,  R.  F.  Harrington,  T.  K.  Sarkar,  and  M.  B.  Bazdar.  This 
* 

program  uses  the  output  of  the  previous  program  as  input .  It  computes 
the  response  of  a  lossless  multiconductor  transmission  line  with  linear 
or  nonlinear  loads.  It  uses  the  modal  time-domain  solution,  which  is  the 
most  accurate  and  f vastest  one.  The  terminating  networks  may  contain 
generators,  resistors,  short  circuits,  open  circuits,  interconnectors, 
inductors,  capacitors,  and  nonlinear  resistors.  The  output  consists  of 
voltage  waveforms  on  all  lines  at  specified  points  along  the  line. 

3)  "Scattering  Parameters  of  Microwave  Networks  with  Multiconductor 
Transmission  Lines,"  by  A.  R.  D jordjevic,  M.  B.  Bazdar,  G.  M.  Vitosevic, 
R.  F.  Harrington,  and  T.  K.  Sarkar.  This  program  uses  the  output  of  the 
first  program  as  input.  It  analyzes  networks  consisting  of  multiconduc¬ 
tor  transmission  lines  and  discrete  elements.  The  output  is  the  scatter¬ 
ing  matrix,  the  impedance  matrix,  and  the  admittance  matrix  of  the 


network.  The  transmission  lines  are  treated  by  modal  analysis  in  the 
time  domain/  including  losses.  The  output  is  in  a  format  similar  to  that 
obtained  from  network  analyzers  used  for  experimental  measurements. 

III.  PLANS  FOR  FUTURE  WORK 

A)  The  work  on  mode  dispersion  on  a  microstrip  transmission  line, 
item  IIB  above,  will  be  published. 

B)  Study  of  the  solution  to  three-dimensional  printed  circuits,  item 
IIC  above,  will  continue.  It  is  planned  to  develop  a  user-friendly 
computer  program  to  handle  three-dimensional  integrated  circuits  of 
arbitrary  shape. 

C)  It  time  permits,  the  full-wave  solution  of  a  microstrip 
transmission  line,  item  IIB  above,  will  be  generalized  to  a 
multiconductor  line  with  conductors  of  arbitrary  shape. 

i 
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Abstract  — The  equivalent  circuit  of  a  via  which  connects  two  semi- 
inf  initvly  long  transmission  lines  through  a  circular  hole  in  a  ground  plane 
is  considered.  The  v-type  equivalent  circuit  consists  of  two  excess  capaci¬ 
tances  and  an  excess  inductance.  They  are  quasi-static  quantities  and  thus 
are  computed  statically  by  the  method  «!  moments  from  the  integral 
equations.  The  integral  equations  are  established  by  introducing  a  sheet  of 
magnetic  current  in  the  electrostatic  case  and  a  layer  of  magnetic  charge  in 
the  magnetostatic  case.  Parametric  plots  of  the  excess  capacitances,  the 
excess  inductance,  and  the  characteristic  admittance  of  the  via  are  given 
for  reference. 

I.  Introduction 

THE  GEOMETRY  of  the  problem  to  be  considered  in 
this  paper  is  shown  in  Fig.  1.  Two  semi-infinitely  long 
transmission  lines,  wire  1  and  wire  2,  are  connected  by  a 
via  through  a  hole  in  a  conducting  ground  plane.  The  via 
consists  of  wire  3  and  wire  4.  The  radii  of  wires  1,  2,  3.  4, 
and  the  hole,  denoted  a,,  a2,  a3,  a},  and  a5,  respectively, 
are  very  small  compared  to  the  heights,  h{  and  Jt2,  of  wire 
1  and  wire  2  with  respect  to  ground.  The  media  in  the 
upper  region  (region  a)  and  the  lower  region  (region  b) 
may  be  different.  Let  us  assume  that  (€t,  /a1  )  and  ( « 2.  ) 

are  the  constitutive  constants  for  region  a  and  region  b, 
respectively.  Also,  the  media  and  all  the  conductors  are 
perfect  (lossless).  For  simplicity,  the  equivalent  circuit  of 
the  via  is  assumed  to  be  w-type,  as  is  shown  in  Fig.  2.  In 
Fig.  2,  Yn |  is  the  characteristic  admittance  of  wire  1  above 
the  ground  plane  and  *02  is  the  characteristic  admittance 
of  wire  2  below  the  ground  plane.  The  circuit  of  Fig.  2  is 
valid  when  only  a  small  portion  of  the  line  voltage  is 
dropped  across  Le  and  when  only  a  small  portion  of  the 
line  current  is  shunted  through  CrX  and  Ct2.  We  desire  to 
determine  the  capacitances  Crl  and  Ct2  and  the  inductance 
Lr.  The  problem  described  here  is  of  practical  interest.  For 
example,  printed  circuits  on  different  sides  of  a  ground 
plane  inside  computers  are  often  connected  by  a  via 
through  a  hole  in  the  ground  plane.  The  related  problems 
of  the  connection  of  two  perpendicular  strips  above  a 
ground  plane  (1),  the  connection  of  two  parallel  wires 
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Fig.  2.  The  equivalent  circuit  for  the  problem. 


above  a  ground  plane  (2],  and  the  connection  of  two 
parallel  strips  above  a  ground  plane  (3),  [4]  were  previously 
considered. 

The  equivalent  capacitance  Ca  of  the  portion  of  the  via 
above  the  ground  plane  is  a  quasi-electrostatic  quantity 
and  is  defined  as  [2],  [3] 


Q  =  lim 

I,  —  00 


Qi +  Q)  A?oi 
v 


Here  Qx  is  the  total  electric  charge  on  the  portion  of  wire  1 
of  length  /,,  Q 3  is  the  total  charge  on  wire  3,  V  is  the 
constant  voltage  maintained  at  the  surface  of  the  wires 
with  respect  to  ground,  and  <?„,  is  the  uniform  charge 
density  on  wire  l  far  away  from  the  via,  or  equivalently, 
the  charge  density  required  to  raise  the  potential  of  wire  1 
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lO'V  volts  if  the  hole  were  closed,  wire  3  were  removed, 
and  wire  1  were  extended  to  infinity.  Henceforth  in  the 
electrostatic  case,  charge  density  means  charge  per  unit 
length.  However,  since  the  numerator  of  the  right-hand 
side  of  (1)  approaches  a  constant  that  is  very  small  com¬ 
pared  to  Q{  and  /t^0l,  as  /,  becomes  large,  numerical 
calculation  of  from  (1)  would  result  in  significant 
error.  To  avoid  this,  we  subtract  the  uniform  charge  den¬ 
sity  from  the  total  charge  density  so  that  the  difference, 
called  the  excess  charge  density,  is  the  unknown  in  the 
boundary  integral  equations.  The  equivalent  capacitance 
Cfl  is  then  the  sum  of  the  total  excess  charge  on  wire  1  and 
wire  3  if  the  voltage  V  is  set  to  one  volt.  Hence  Crl  is  also 
called  the  excess  capacitance  of  the  upper  part  of  the  via 
(wire  3).  In  a  similar  manner,  the  equivalent  capacitance 
Ct2  >s  defined  and  is  called  the  excess  capacitance  of  the 
lower  part  of  the  via  (wire  4). 

The  equivalent  inductance  Le  is  a  quasi-magnetostatic 
quantity  defined  by  [2],  [4] 
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Fig  3.  The  electrostatic  problem  divided  into  uu>  parts  (u|  llicelcctiic 
field  remains  unchanged  in  region  a  if  the  hole  is  closed  and  M  is 
placed  above  the  hole,  (b)  The  electric  field  remains  unchanged  in 
region  b  if  the  hole  is  dosed  and  -  M  is  placed  just  below  the  hole 

and  region  b  if 

34 >*  S4>* 

e .  — —  =  e ,  in  the  hole  (4) 

(?/»  '  thi 

where  <#>'  is  the  electric  potential  in  region  /,  /  =  a,b.  In 
addition  to  (4)  the  following  boundary  conditions  must 
also  be  satisfied: 


where  Ai  is  the  total  magnetic  vector  potential  on  wire  i. 
/  =-1,2,3, 4,  due  to  the  steady  electric  current  of  filamen¬ 
tary  strength  /  flowing  from  wire  1  through  the  via  to  wire 
2.  ,410  is  the  uniform  magnetic  vector  potential  on  wire  1 
far  away  from  the  via  or,  equivalently,  the  magnetic  vector 
potential  that  would  be  produced  by  an  /  ampere  current 
on  wire  1  if  the  hole  were  closed,  wire  3  were  removed,  and 
wire  1  were  extended  to  infinity.  A  20  is  similarly  defined. 
The  sum  of  the  four  integrals  in  (2)  is  a  line  integral  from  a 
point  far  to  the  left  on  the  surface  of  wire  l  to  a  point  far 
to  the  right  on  the  surface  of  wire  2.  Equation  (2)  is 
derivable  from  [4,  eq.  (1)J.  Le  is  also  called  the  excess 
inductance  of  the  via.  For  convenience,  we  call  the  first 
bracketed  term  the  excess  inductance  of  the  upper  via 
(wire  3)  and  the  second  bracketed  term  the  excess  induc¬ 
tance  of  the  lower  via  (wire  4). 

II.  Formulation 

The  excess  capacitances  Q  and  C, 2  are  quasi-electro- 
static  quantities  and  thus  computed  in  the  electrostatic 
case.  As  shown  in  Fig.  3,  we  first  close  the  hole  by  a 
conductor  and  place  a  sheet  of  magnetic  current  M  just 
above  the  hole  and  -  M  just  below  it.  Steady  in  that  it  has 
no  surface  divergence,  this  magnetic  current  is  related  to 
the  electric  field  EA  over  the  hole  by 

M~EAXn  (3) 

where  n  is  the  unit  vector  normal  to  the  ground  plane  and 
pointing  upwards  (from  region  b  to  region  a).  By  the 
uniqueness  theorem  (3,  sec.  3-3]  and  the  equivalence  prin¬ 
ciple  (S,  sec.  3-3],  the  field  remains  the  same  in  region  a 


4>a  =  V  on  wires  1  and  3 

4>b  =  V  on  wires  2  and  4. 


(5) 


Neither  wire  3  nor  wire  4  is  connected  to  the  ground  plane 
in  Fig.  3. 

The  total  charge  density  may  be  recognized  as  the  sum 
of  the  excess  charge  density  and  uniform  charge  density, 
i.e., 

I  <7,i  +  <7m  on  1 

I  <7,3  on  wire  3 

q  =  {  .  (6) 

\  qrt  on  wire  4 

W,2  +  <7»2  on  wire  2 


where  the  subscript  e  on  q  denotes  excess  charge  density. 
(701  and  ql)2  are  known  and  are  given  by 

2w<,F 

q"x=  in(2A,/fl,)’ 

2m2V 

q"2=  \n{2h2/a2)' 

If  the  potential  V  is  set  to  one  volt,  the  total  excess  charge 
will  give  rise  to  the  desired  excess  capacitance,  i.e.. 


/i ,  »  a , 


/',  »  a  2. 


(7) 


Q-/  <7,i  d/  +  f  q,yJI 

•'wire  1  •'wire  3 

C,2  =  /  <7,2  <7/ +  /  q,<dl.  (8) 

•'wire  2  ■'wire  4 

Let  4>(qe, )  denote  the  potential  due  to  </,,,  in  the  pres¬ 
ence  of  the  completed  ground  plane  (hole  closed),  4>[  (M) 
the  potential  in  region  a  due  to  M  residing  on  the  region  a 
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side  of  the  completed  ground  plane,  and  <p  ;  M )  the 
potential  in  region  b  due  to  M  residing  on  the  region  h 
side.  Then, 


=»  <p(qm )  +  <f>(</»i )  +  )  +  4>  ‘  (  M  ) 

<f»h  =  <f> (</<«)  +  ^  <M</.  ^  (M).  (9) 

Note  that  the  electric  potential  due  to  a  sheet  of  steady 
magnetic  current  is  analogous  to  the  magnetic  scalar 
potential  due  to  a  sheet  of  steady  electric  current.  The 
latter  is  studied  in  many  fundamental  electromagnetic  field 
theory  books,  for  example,  (6).  Substitution  of  (9)  into  (4) 
and  (5),  with  V  being  set  to  one  volt,  yields  integral 
equations  for  the  excess  charge  densities: 

$(<7»i )  +  <#>(  <7,j )  +  <f» f  ( M ) 

=  1 -<#>(</„,)  on  wires  1  and  3 

$(<7,2) +  ^(<7,4)“  (M) 

“  1  -  </> ( t/o2 )  on  wires  2  and  4 
H(qti)  d4>(qe  3)  d<t>*(M) 

<)<t>(qe 4)  d<>  ( M ) 

2  dn  2  dn  <2  dn 


-“«i 


d<t>(qoi) 

dn 


+  «2 


d<ft(<7o2) 

dn 


in  the  hole.  (10) 


Although  qel  and  qe 2  exist  on  the  semi-infinitely  long 
wires  1  and  2.  they  decay  to  zero  rapidly  as  one  moves 
away  from  the  via.  We  may  truncate  qri  and  qe,  and  the 
boundary  equations  on  wire  1  and  wire  2  at  some  dis¬ 
tances,  say  3 hl  and  3 h2,  from  the  via.  By  doing  so,  we 
neglect  the  contribution  to  the  excess  capacitances  from 
the  excess  charge  beyond  the  lengths  37i,  on  wire  1  and  3h2 
on  wire  2.  Now,  the  method  of  moments  may  be  used  to 
solve  for  qrl,  qe2,  qeJ,  qr4,  and  M  numerically.  We  divide 
wire  1  and  wire  2  (truncated)  and  wire  3  and  wire  4  into 
subsections,  assume  uniform  charge  distribution  on  each 
subsection,  and  enforce  the  first  two  of  equations  ( 10)  at 
the  center  of  each  subsection.  Furthermore,  we  divide  the 
hole  into  annuluses,  assume  uniform  circulating  current 
distribution  on  each  annulus,  and  enforce  the  third  of 
equations  (10),  which  is  averaged  over  the  interval  from  0 
to  2 it  for  the  azimuthal  variation,  at  the  midpoint  between 
the  edges  of  each  annulus.  A  detailed  discussion  of  the 
moment  method  as  applied  to  this  problem  is  presented  in 

17]. 

Now,  we  turn  to  the  magnetostatic  case  to  compute  the 
excess  inductance.  Similar  to  the  electrostatic  case,  we 
close  the  hole  by  a  conductor  and  place  a  layer  of  magnetic 
charge  density  m  just  above  the  hole  and  -  m  just  below 
it,  as  is  shown  in  Fig.  4.  This  m  is  equal  to  the  normal 
magnetic  field  over  the  hole.  Again  by  the  uniqueness 
theorem  and  the  equivalence  principle,  the  magnetic  field 
remains  unchanged  in  region  a  and  region  b  if 

(ID 

is  enforced,  where  /#,'an  is  the  tangential  magnetic  intensity 


1  hi 


I'ig.  4  The  magnclovlalic  problem  divided  into  two  parts,  (a)  The 
magnetic  field  remains  unchanged  in  region  a  if  the  hole  is  closed  and 
m  is  placed  just  above  the  hole,  (b)  The  magnetic  field  remains 
unchanged  in  region  b  if  the  hole  is  closed  and  -  ni  is  placed  just 
below  the  hole. 

over  the  hole  in  region  i,  i  =  a.b.  Let  //  (wire  i)  denote 
the  magnetic  intensity  due  to  the  electric  current  on  wire  / 
in  the  presence  of  the  completed  ground  plane  for  /  = 
1.2.3, 4,  //'  («t)  the  magnetic  intensity  due  to  m  above 
the  completed  ground  plane,  and  W“(m)  the  magnetic 
intensity  due  to  m  below  the  completed  ground  plane.  We 
can  write 

H*  -  //(wire  1)+  //(wire  3)+  H *  (m) 

//h  =  //(wire  2)  +  //(wire  4)  -  H~  (m ).  (12) 

Note  that  the  magnetic  intensity  due  to  a  layer  of  magnetic 
charge  is  analogous  to  the  electric  field  due  to  a  layer  of 
electric  charge,  as  discussed  in  [6|.  Therefore,  H*  (m)  and 
//  (m)  may  be  represented  by  the  gradients  of  some 
scalar  functions  vk  ’  ( «i )  and  'F  '  ( m ).  That  is. 

//  *  ( m )  =  -  V  * 4  ( m ) 

//  (m)  =  — V*  (m).  (13) 

Substitution  of  (12)  and  (13)  into  (II)  gives 

V'F,(w)  +  V'F  (/»)|,an 

=  //(wire  1 )  -  //(wire  2)  +  //(wire  3)  -  //(wire  4)|lin. 

(14) 

Equation  (14),  an  integral  equation  for  m,  implies  that 
'F  ’  (m)+  ♦  ‘  (m)  equals  'F"K.  where  'F"K  is  a  potential 
whose  gradient  is  the  right-hand  side  of  (14).  To  com¬ 
pensate  for  the  fact  that  'Finc  is  only  known  to  within  an 
additive  constant,  we  require  that  the  total  magnetic  charge 
associated  with  m  vanish.  The  moment  method  may  be 
applied  now  to  solve  the  scalar  equation  derived  from  (14) 
subject  to  the  above  constraint  on  m. 

Since  there  is  no  coupling  between  wire  1  and  wire  3  and 
between  wire  2  and  wire  4,  we  have 

/4,  =  /4,(wire  1)+  At'  ( m ) 

A 2  =  ,42(wire  2)-  A{  (m) 

A ,  =  A ,  ( wire  3)  +  A  ( m ) 

A4  ~  A 4 (wire  4)-  A4  (mi) 


(15) 
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where  i4,(wire  i)  is  the  magnetic  vector  potential  on  the 
surface  of  wire  /  due  to  the  current  /  on  wire  i  in  the 
presence  of  the  completed  ground  plane.  /  =  1,2.3,4. 
Moreover,  A* (m)  is  the  vector  potential  on  the  surface  of 
wire  f,  /  —  1 , 3,  due  to  m  residing  on  the  region  a  side  of 
the  completed  ground  plane,  and  A,  ( m )  is  the  vector 
potential  on  the  surface  of  wire  i,  ;  =  2.4.  due  to  m 
residing  on  the  region  b  side.  Note  that  (15)  is  valid  only 
for  the  vector  potential  component  tangent  to  each  wire 
and  that  the  other  component  is  not  of  interest.  Letting  / 
be  one  ampere  and  putting  (15)  into  (2),  we  get  the  excess 
inductance  of  wire  3: 

Lel  *  J  (/4,(wire  1)- Al0)-dl  +  f  A  j(wire  3)  dl 

•'wire  1  •'wire  3 

+  [  Ai(m)  dl+(  Aj  (m)- dl.  (16) 

*  wire  1  ■'wire  3 

The  excess  inductance  of  wire  4  is  obtained  from  the  above 
expression  by  replacing  all  the  1  ’s  by  2’s,  3’s  by  4's,  m ’s  by 
—  m ’s,  and  superscripts  +  by  - ,  that  is. 

Le2  **  f  (<42(wire  2)- Am)dl  +  (  /l4(wire  4) •  dl 

•'wire  2  •'wire  4 

-[  A2(m)dl-f  A;(m)dl.  (17) 

Because  of  the  small  radii  of  the  wires,  the  current  on  the 
surface  of  each  wire  may  be  approximated  by  a  filamen¬ 
tary  current  /  on  the  axis  of  the  wire.  Thus  the  first  two 
integrals  in  (16)  can  be  evaluated  analytically  (2|.  The  last 
two  integrals  may  be  expressed  in  terms  of  magnetic 
intensities  according  to  Stokes's  theorem.  The  result  is 


Zh,  In 


~4h,  +  ut  +  a  3 


■/ «:( 

•V, 


m)  dt  (18) 


in  the  radial  direction  with  the  variation  being  of  the  form 


p\j\  ~(p/a<)~ 

where  p  is  the  radial  distance  from  the  center  of  the  hole 
and  K  is  a  constant  determined  by 


/ dl  =  1 . 


The  integration  path  is  chosen  in  the  radial  direction  from 
the  surface  of  wire  3  to  the  edge  of  the  hole.  Substitution 
of  (19)  into  (20)  gives 


ln(  +  /(05  ’j)2-  1  ) 


where  in  the  last  integral,  5,  is  the  planar  surface  bounded 
by  the  axes  of  wire  1  and  wire  3  and  the  ground.  The 
subscript  n  on  H*(m)  denotes  the  component  normal  to 
S,  and  pointing  into  the  paper.  The  expression  for  l.r2  is 
similar  to  (18)  in  form.  The  excess  inductance  of  the  via  is 
then  the  sum  of  Ltl  and  L,2. 


III.  Numerical  Results  and  Discussion 

In  order  to  implement  the  computer  program  on  a 
PC/AT,  we  simplify  the  computation  by  choosing  only 
one  expansion  function  for  the  magnetic  current  and  none 
for  the  magnetic  charge.  This  is  justified  as  follows.  Since 
the  radius  of  the  hole  is  very  small  compared  to  the  heights 
A,  and  A2,  the  couplings  (electric  and  magnetic)  through 
the  hole  between  wire  1  and  wire  2  are  negligible  to  the 
first-order  approximation.  In  other  words,  the  fields  in  the 
hole  are  primarily  from  the  via  (wire  3  and  wire  4).  In  the 
electrostatic  case,  the  tangential  electric  field  in  the  hole  is 


Thus  the  magnetic  current  is  circulating  and  its  amplitude 
is  given  by 

1  1 

Af - 7 - - r--r— — J-  (22) 

In  |  <is/ri,  + /(«,/«,)  -Ij  pyl  ~(p/<i,) 

In  the  magnetostatic  case,  the  normal  magnetic  field  in 
the  hole  is  negligible.  Hence  the  magnetic  charge  may  be 
neglected  and  a  closed  form  for  the  excess  inductance  of 
the  via  is  obtained.  That  is, 

<251 

where  x,  and  k2  are  constants.  Approximately,  k,  =  k2  = 
0.5413. 

In  Figs.  5-7.  the  normalized  excess  capacitance  of  the 
upper  via  (wire  3)  C(.|/(t|<J3)  is  plotted.  The  curves  are 
also  applicable  to  the  normalized  excess  capacitance 
C,2/(c2a,)  of  the  lower  via  (wire  4)  if  all  the  subscripts  1 
are  replaced  by  2.  In  Figs.  8  and  91  the  normalized  excess 
inductance  Lr/(2ju,«,)  of  a  symmetric  via  is  plotted.  A  via 
is  symmetric  if  (It,.  p,)  =  (h2,  a2,e2,  p2).  These  plots 

can  also  be  viewed  as  plots  for  the  normalized  excess 
inductance  /(p,n,)  of  the  upper  via.  If  the  subscripts  1 
arc  replaced  by  2.  they  become  plots  of  l.r2/(ii2<i , ),  the 
normalized  inductance  of  the  lower  via.  In  Figs.  10  and  1 1 
the  characteristic  admittance  r),K,  of  a  symmetric  via  is 
plotted.  Here  tj,  is  the  intrinsic  impedance  of  the  medium. 
The  characteristic  admittance  of  the  via  is  defined  as 


r,  =  {(Crl  +  Cr2)/Lf.  (24) 

If  a  TEM  wave  approaches  the  via  along  wire  1  and  if 
wire  2  is  terminated  with  a  matched  load  ( >'02 ).  then  it  can 
be  shown  that  the  least  reflection  will  occur  when 

Y,  ~  •  (25> 

*  Note  that  the  curses  in  Pigs.  X  and  9  arc  plotted  under  the  assumption 
that  the  hole  is  small  in  relation  to  the  heights  /i,  and  li,.  To  emphasize 
this  and  to  he  consistent  » till  other  plots,  we  add  the  restriction  u,  2 n, 


«  - i  m  .  in - »»  - jn - >  » 


Fig.  6.  Normalized  excess  capacitance  of  the  upper  via  (<Jt  -a  j.  k  -  Fig.  8.  Normalized  excess  inductance  of  a  symmetric  via  (a5  -  2a5, 

“s  /<*})■ 

To  see  this,  wc  notice  that  the  voltage  reflection  coefficient  at  the  T  plane  on  wire  1  in  Fig.  2  is  given  by 


.  i 

j  w2f  rC,2 

l  2  J 

J»L.  "“‘t. 

Wl 

’  K01  ( 1  -  w2L,C,2)  +  K02(l  -  W2L,C„) .  1 

r 

L  2} 

j»L.  +  r°'Y°1  + 

l  2 

Ct2y  2  ) 

where  a  is  the  angular  frequency.  Usually,  u2Lt{Ctl  + 
C,2)  <  1.  The  above  expression  then  reduces  to 

„  rn  ~  *02  +  juLe  Y0l  Y02  -  ju  ( C„  +  C,i ) 


To  design  a  reflectionless  via,  consider  the  following 
numerical  example.  Suppose  that 


p  roi  ~  *02  +  JuLt roi ro2~  +  t-'j)  Ai  =  A2  =  1.00cm 

*ot +  *02  +  y«f»»* oi*o2 +  JW(C*  i +  C2)  a,  =  a2  =»  0.10  cm 

When  (25)  is  satisfied,  the  magnitude  of  the  numerator  of  <j5  =  2a 3  (28; 

(27)  assumes  its  minimum,  the  magnitude  of  the  de¬ 
nominator  is  roughly  T01  +  Y0J,  and  thus  IT)  is  minimized.  and  that  we  wish  minimize  reflection  from  the  via.  Since 
The  reflection  from  the  via  is  minimized.  Furthermore,  if  2ff  2  10 

(25)  is  satisfied  and  if  the  system  is  symmetric  ( K01  *  *02)*  *01  =  *02=  '  7  771— :  7-—  (29 

there  is  no  reflection  from  the  via.  That  is,  all  of  the  power  ijln(2At/ ax)  tj 

from  the  incident  wave  will  be  transmitted  through  the  via  the  refiectionless  condition  becomes 
to  the  matched  load.  In  this  case,  the  via  is  called  reflec- 

tionless.  ^,  =  2.10.  (30 
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Fia.  9.  Normalized  excess  inductance  of  a  symmetric  via  (a,  -  2uj. 
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Fig.  10.  Normalized  characteristic  admittance  of  a  symmetric  via  (a, 
2aj,  *  -a, /a,). 
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Fig.  11.  Normalized  characteristic  admittance  of  a  symmetric  via  (a}  — 
2aJt  k  «■  6,/a,). 

Using  the  ft  -10  curve  in  Fig.  11,  we  find  that  the  above 
condition  is  satisfied  when  a,  —  0.034  cm.  Hence,  given 
[28),  the  via  will  be  reflectioniess  when  the  radius  of  the 
ria  is  0.034  cm. 
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The  Inductance  Matrix  of  a  Multiconductor 
Transmission  Line  in  Multiple  Magnetic  Media 

JOSEPH  R.  MAUTZ,  senior  member,  ieee, 

ROGER  F.  HARRINGTON,  fellow,  ieee,  and 
CHUNG-I G.  HSU,  student  member,  ieee 

Abstract  —Co mUh  a  nuMconductor  transmission  lino  consisting  ol  Ne 
conducting  cylinders  in  inhomogeneous  media  consisting  ol  N ,  homoge¬ 
neous  regions  with  permeabilities  p,  and  permittirities  The  inductance 
matrix  ( L]  for  the  line  is  obtsined  by  solving  the  magnetostatic  problem  of 
Vf  conductors  in  Nj  regions  with  permeabilities  p,.  The  capecitance 
matrix  (C|  for  the  line  is  obtsined  by  solving  the  eiectrostatic  problem  of 
conductors  in  Sj  regions  with  permittirities  It  is  shown  that 
[Z-l -PoToIC']-1.  where  (C'|  is  the  capacitance  matrix  of  an  auxiliary 
electrostatic  problem  of  Nc  conductors  in  Nj  regions  with  relative  permit¬ 
tivities  set  equal  to  the  reciprocals  of  the  relative  permeabilities  of  the 
magnetostatic  problem,  i.e.,  <;/<0  -  p0/p,. 

I.  Introduction 

Fig.  1  shows  a  multiconductor  transmission  line  consisting  of 
iV  conductors  and  Nlt  insulating  materials  above  a  perfectly 
conducting  ground  plane.  The  system  is  uniform  in  the  z  direc¬ 
tion  (direction  of  propagation),  and  the  cross-sectional  shapes  of 
the  conductors  and  insulators  are  arbitrary.  The  insulators  have 
arbitrary  permeabilities  p,  and  permittivities  An  upper  ground 
plane  could  be  present  and  treated  as  an  additional  conductor  in 
a  manner  similar  to  that  of  [1]. 

To  the  quasi-static  approximation,  the  multiconductor  trans¬ 
mission  line  is  characterized  by  a  capacitance  matrix  [C],  ob¬ 
tained  from  an  electrostatic  analysis,  and  an  inductance  matrix 
[  L],  obtained  from  a  magnetostatic  analysis.  The  formulation  of 
the  problem  for  [CJ  and  a  numerical  algorithm  for  its  computa¬ 
tion  are  given  in  (1).  In  ihe  appendix  of  [1]  it  is  shown  that  if  the 
insulating  matter  is  nonmagnetic,  the  inductance  matrix  of  the 
line  is  given  by 

[L]“Po*o[^ol  (1) 

Here  [Q]  is  the  capacitance  matrix  of  the  multiconductor  trans¬ 
mission  line  if  all  dielectric  constants  are  set  equal  to  1. 

In  this  paper  we  shall  show  that  when  the  insulators  are 
magnetic,  a  modified  relationship  holds.  In  particular,  it  is 

m-^oicr'  (2) 

where  [C]  is  the  capacitance  matrix  of  the  line  if  all  relative 
permittivities  are  set  equal  to  the  reciprocals  of  the  relative 
permeabilities,  i.e., 

— .  « -1,2,  •••,*„.  (3) 

*0  A, 

Note  that  the  auxiliary  problem  for  Finding  [C']  has  relative 
permittivities  less  than  unity  if  the  corresponding  relative  per¬ 
meabilities  are  greater  than  unity. 

Relationship  (2)  is  strictly  true  only  if  all  conductors  are 
perfect  It  assumes  that  all  current  in  the  magnetostatic  problem 
flows  on  the  surfaces  of  conductors.  For  actual  conductors, 
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Fig.  1.  A  multiconductor  transmission  line  in  a  multilayered  dielectric  and 
magnetic  region  above  a  ground  plane. 


depending  on  the  frequency,  some  of  the  current  flows  internal  to 
the  conductors.  Relationship  (2)  is  then  only  approximate. 

Although  (2)  can  be  inferred  from  the  last  two  sentences  in 
[2,  sec.  V],  we  have  not  seen  an  explicit  proof  in  the  literature. 
The  purpose  of  this  short  paper  is  to  give  a  simple  proof  of  (2) 
for  the  multiconductor  case  with  layered  media. 

II.  Proof  of  (2) 

e- 

The  electrostatic  problem  from  which  [C]  is  calculated  is 
formulated  in  detail  in  [1].  We  shall  refer  to  that  formulation 
when  needed.  The  magnetostatic  problem  from  which  [L]  is 
calculated  is  an  extension  of  the  formulation  for  nonmagnetic 
media,  given  in  the  appendix  of'[l].  The  formulation  for  magnetic 
media  is  given  below. 

The  inductance  matrix  [L]  is  an  Nc  X  Sc  matrix  that  satisfies 

+  -[L\1  (4) 

where  $  and  Tare  V,  x  1  column  vectors.  The  yth  element  of  f 
is  the  z-directed  conduction  current  on  the  jth  conductor.  The 
/  th  element  of  \i>  is  the  x-directed  magnetic  flux  passing  between 
a  unit  length  of  the  i  th  conductor  and  the  lower  ground  plane.  If 

'■rr(c'i*  (5> 

Po*o 

and  if  [C']_1  exists,  then  (4)  will  imply  that  the  desired  relation¬ 
ship  (2)  is  true.  In  the  remainder  of  this  section,  we  establish  (5). 

We  formulate  the  magnetostatic  problem  with  magnetic  media 
in  terms  of  the  total  electric  current  Jru.  on  the  surfaces  of  the 
conducting  cylinders  and  on  the  boundaries  between  different 
magnetic  media.  The  total  electric  current  on  the  conducting 
surfaces  is  the  conduction  current  plus  the  magnetization  current 
The  total  electric  current  on  the  magnetic  media  boundaries  is 
the  magnetization  current  The  magnetic  flux  density  B  is  given 
by 

B-v  XA-VAtXM.  (6) 

where  A.  is  the  only  component  of  A  due  to  steady  current 
flowing  in  the  z  direction.  For  the  two-dimensional  problem,  it  is 
given  by 

(7) 

where  /,  denotes  the  j  th  interface.  The  first  Nc  interfaces  are  the 
surfaces  of  the  Nc  conductors.  If  the  upper  ground  plane  is 
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present,  the  (Se  +  l)th  interface  is  the  surface  of  this  plane.  The 
Inst  NJ  jn^rfy**  are  the  magnetic  media  boundaries.  Thus, 

M-Mx  +  Nj  (8) 

when  is  Nt  it  the  upper  (round  plane  is  absent  If  the  upper 
ground  plane  is  present,  then  Mx  is  Nc  + 1.  The  number  Nj  of 
magnetic  media  boundaries  could  be  greater  than  - 1  because 
these  boundaries  are  arbitrarily  shaped,  not  necessarily  parallel  to 
one  another  or  to  the  lower  ground  plane.  In  (7),  p'  is  the 
position  vector  of  dl’,  and  p'  is  the  position  vector  of  the  image 
of  dl'  about  the  lower  ground  plane. 

Substituting  (7)  into  (6)  and  a«»ming  that  p  is  not  on  any  of 
the  interfaces  (/, },  we  obtain 


>'\2  Ip  -  pT 


xu.dl'.  (9) 


The  limits  of  (9)  as  p  approaches  /,  from  either  side  are 


XU,  dl' 


Pojt(p) 


Ip  on  I, 

\  /  “1,2,-  •  i 


where  f,  denotes  the  principal  value  of  the  integral  over  and 
m  is  a  unit  vector  normal  to  I,  at  p.  Moreover,  B*  (p)  is  B(p)  on 
the  side  of  I,  toward  which  *  points,  and  B~  (p)  is  B( p)  on  the 
other  side  of  /,. 

On  the  surface  of  the  ith  conductor.  A,  is  constant  and  is  the 
x-directed  magnetic  flux  passing  between  a  unit  length  of  the 
<  th  conductor  and  t>.e  lower  ground  plane.  Hence,  similar  to 
[I,  eq.  (A2)l, 

^(p)~i,<  l/S ?t...  m  (ll> 


If  the  i  th  conductor  is  of  finite  cross  section,  then  /,  is  a  closed 
curve  on  which 

,,;r(p).-»xi*(p)  (15) 

Po 

•Mp  )--r7-7»x**(p)  (16) 

p,  (p) 

where  it  is  the  unit  normal  vector  that  points  outward  from  the 
surface  of  the  conductor.  Furthermore,  B*  (p)  and  p*  are, 
respectively,  the  magnetic  flux  density  and  the  permeability  just 
outside  the  conductor.  In  (16),  Jc(p)  is  the  conduction  current 
on  the  conductor.  Equations  (IS)  and  (16)  imply  that 

<17) 

on  the  surface  of  the  i  th  conductor,  provided  this  conductor  is  of 
finite  cross  section. 

If  the  i  th  conductor  is  an  infinitesimally  thin  strip,  then  l,  runs 
from  one  edge  of  the  strip  to  the  other  on  which 

<U) 

Substitution  of  (10)  for  B1  (p)  in  (18)  leads  to 


,  ,  .  Po 


P,*(P)  +  P.“ 


Mp) 


If  the  upper  ground  plane  is  present,  it  is  the  (Aff  +  l)th  conduc¬ 
tor  and,  because  A,  vanishes  on  it,  ^  t  -  0.  Substitution  of  (7) 
into  (11)  gives 

Po  £  r  .  ,  (Ip~pI\  ...  .  (P°a/. 

2^?iV(P)  ll  P-P'lJ  '■ 

(12) 

Continuity  of  the  tangential  component  of  magnetic  intensity 
on  the  magnetic  media  boundaries  requires  that 

B*(  p)  B"(p)1  [pon/j 

"X[  P*  ~  h  ]  ’ 

(13) 

where  n,  B*( p),  and  B~(p)  are  the  same  as  in  (10).  Further¬ 
more,  p*  is  the  permeability  on  the  side  of  l,  toward  which  n 
points  and  is  the  permeability  on  the  other  side  of  /,. 
Substituting  (10)  into  (13)  and  then  dividing  (13)  by  u,[l/p*  - 
1/pT  ],  we  obtain 

P±_ _ fb]  .  /  _v  ,  Po  £  t  ,  / 


2w[p;(p)  m,*(p) 

(  p  ~ p  p  -  p'  \ 
\|p-p'l2  Ip  —  P'l2  / 


E  fMp') 

j-l  'i 


Mo  M/  Mi  f  /  \  ,  Mo  p  f »  /  /\ 

t  T7T  /t(,) 

. J 

•f  p~.~  |  ndl’  -0, 

l  IP  -  P  i  IP  -  P I  / 

I  pool, 

\i-A/i  +  l,Af,+2.--,M.  (14) 


IIP-P'I2  IP -Pi2/ 

on  the  surface  of  the  i  th  conductor,  provided  this  conductor  is  of 
zero  thickness. 

Now,  consider  the  auxiliary  electrostatic  problem  which  has 
the  same  geometry  as  that  of  the  present  magnetostatic  problem, 
but  with  relative  permittivities  e,'/«0  s®1  equal  to  Po/Pr  The 
formulation  presented  in  (1)  is,  in  fact,  valid  for  dielectric  media 
of  arbitrary  shape.  The  unit  vector  uv  in  (11)  of  [1]  should  be 
replaced  by  »  when  the  dielectric  media  arc  arbitrarily  shaped.  It 
is  clear  that  (12),  (14),  (17),  and  (19)  have  the  same  mathematical 
forms  as  (9),  (11),  (15),  and  (17)  of  [1J,  respectively.  Therefore, 
the  solution  of  the  magnetostatic  problem  can  be  related  to  that 
of  the  auxiliary  electrostatic  problem  by 

4°{P ') - - •Hp').  (20) 

where  J£'\p’)  is  the  conduction  current  of  the  magnetostatic 
problem  when  -1  is  the  only  nonzero  magnetic  flux  and 
oJ.'Hp')  is  the  free  charge  of  the  auxiliary  electrostatic  problem 
when  the  potential  of  the  r  th  conductor  is  unity  and  all  other 
conductors  are  grounded.  Multiplying  (20)  by  and  summing 
over  i,  we  obtain 

EW)*-—  E«*°(p')*,-  (2i) 

,-l  Po*Oi-l 

After  noting  that  the  right-hand  side  of  (21)  is  7t (p'),  wc  in¬ 
tegrate  (21)  over  l,  to  obtain  (5)  with  the  y/th  element  of  [C 1 
given  by 


c;~jo^(p')  <u\ 


i.j- 1.2, •••,*(. 


(22) 
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.  ‘Premultiplicatiofl  of  (5)  by  (C']~l  yields 

+  -W0[C]-ir.  (23) 

The  inverse  of  (C'J  exists  because  (C'j  is  positive  definite,  which 
can  be  concluded  from  the  fact  that  the  electrostatic  energy 
stored  in  the  system  is  always  greater  than  zero  with  nontrivial 
free  charge  distribution  on  the  conductors.  Comparison  of  (23) 
with  (4)  gives  the  desired  relationship  (2). 

III.  Conclusion 

A  simple  relationship  between  the  inductance  matrix  and  the 
auxiliary  capacitance  matrix  has  been  given.  Thanks  to  this 


relationship,  the  computer  code  given  in  (1)  and  (3)  for  obtaining 
the  capacitance  matrix  of  the  electrostatic  problem  can  be  to 
obtain  the  inductance  matrix  of  the  magnetostatic  problem. 
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Boundary  Integral  Formulations  for  Homogeneous  Material  Bodies 
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Abstract-  There  are  many  boundary  integral  formulations  for  the  problem  of  electro¬ 
magnetic  scattering  from  and  transmission  into  a  homogeneous  material  body.  The  only 
formulations  which  give  a  unique  solution  at  all  frequencies  are  those  which  involve  both 
electric  and  magnetic  equivalent  currents,  and  satisfy  boundary  conditions  on  both  tan¬ 
gential  7  and  tangential  7.  Formulations  which  involve  only  electric  (or  magnetic) 
equivalent  currents,  and  those  which  involve  boundary  conditions  on  only  tangential  E 
(or  tangential  7 )  are  singular  at  frequencies  corresponding  to  the  resonant  frequencies 
of  a  resonator  formed  by  a  perfect  conductor  covering  the  surface  of  the  body  and  filled 
with  the  material  exterior  to  the  body  in  the  original  problem. 

I.  INTRODUCTION 

In  electromagnetic  theory,  a  boundary  integral  equation  is  one  which  involves  the 
integral  of  an  unknown  source  times  a  Green’s  function  over  the  surface  bounding 
a  material  body.  In  this  paper  we  consider  boundary  integral  equations  for  deter¬ 
mining  the  electromagnetic  field  internal  and  external  to  a  homogeneous  material 
body,  when  it  is  excited  by  sources  external  to  the  body.  The  problem  when  the 
sources  are  internal  to  the  body,  or  when  sources  are  both  external  and  internal, 
requires  a  minor  modification  of  the  theory. 

There  are  infinitely  many  boundary  integral  equations  that  can  be  formulated 
for  calculating  the  electromagnetic  field  interned  and  external  to  a  homogeneous 
material  body.  In  this  paper  we  formulate  two  classes  of  such  equations,  which 
we  call  source  formulations  and  field  formulations.  In  the  source  formulation  the 
problem  is  formulated  in  terms  of  unknown  electric  and  magnetic  equivalent  cur¬ 
rents  separately  for  the  internal  and  external  regions,  and  the  boundary  conditions 
are  applied  to  tangential  E  and  tangential  H  on  the  boundary  surface  5 .  This 
leads  to  two  equations  in  four  unknowns,  and  further  restrictions  must  be  made  on 
the  sources.  In  the  field  formulation,  a  single  set  of  equivalent  currents  {J,A1} 
is  used  to  produce  the  field  in  both  the  internal  and  external  regions,  and  the 
boundary  conditions  are  again  applied  to  tangential  E  and  tangential  H  on  5 . 
This  leads  to  four  equations  in  two  unknowns,  and  we  must  make  a  further  choice 
of  how  to  satisfy  the  equations.  Some  choices  of  sources  and  field  equations  lead  to 
operator  equations  that  are  singular  at  certain  “resonant  frequencies”,  and  other 
choices  lead  to  operator  equations  that  have  unique  solutions  at  all  frequencies. 

A  number  of  different  boundary  integral  equations  for  homogeneous  material 
bodies  have  been  used  for  computation  by  various  authors.  For  two-dimensional 


3 


Harrington 


problems,  a  formulation  obtained  by  Muller  (1]  has  been  used  by  Solodukhov  and 
Vasil’ev  [2j  and  by  Morita  [3).  A  formulation  equivalent  to  our  JE-field  method 
has  been  used  by  Wu  and  Tsai  [4],  and  by  Arvas,  Rao,  and  Sarkar  [5|.  For  three- 
dimensional  problems,  a  formulation  obtained  by  many  investigators,  called  the 
PMCHW  formulation  in  [6],  has  been  used  by  Wu  (7),  by  Mautz  and  Harrington 
(6]  and  by  Umashankar,  Taflov,  and  Rao  in  [8).  For  references  to  methods  not 
classified  as  boundary  integral  formulations,  see  [8]. 

II.  DEFINITION  OF  OPERATORS 

The  problem  to  be  considered  is  shown  in  Fig.  1.  The  homogeneous  mate¬ 
rial  body  is  bounded  by  5  and  has  constitutive  parameters  £;  (t  denotes 
internal ) .  The  region  external  to  the  body  is  also  homogeneous  and  has  constitu¬ 
tive  parameters  pe,  e«  (e  denotes  external).  The  impressed  sources  are  denoted 
7*mp,  A/,mp ,  and  the  field  that  they  produce  with  the  body  absent  is  denoted 
jE*mp,  /f*mp .  The  unit  normal  pointing  outward  from  S  is  denoted  n .  The  total 
field  in  the  internal  region  is  denoted  E  ,H  .  The  total  field  in  the  external  region 
is 

E  =  Eimp  +  E“at  (1) 

H  =  /Tmp  +  H‘cai  (2) 

where  the  scattered  field  E>cai ,  H*cai  is  that  produced  by  the  material  body. 
Lossy  bodies  can  be  treated  by  considering  Cj ,  or  Hi  >  or  both,  to  be  complex. 


n 


external  region 

/*.,  c. 


I  =  Tmr  +  T'" 
~B  =  F"'  +  h““ 


Figure  1.  A  material  body  bounded  by  5  in  the  impressed  field  produced 
by  rmp,  ATmp . 

We  use  electric  and  magnetic  surface  currents  radiating  into  homogeneous  me¬ 
dia  of  infinite  extent  as  equivalent  sources.  The  notation  Je,  Me  denotes  equiv¬ 
alent  currents  on  5  radiating  into  media  having  constitutive  parameters  /re,  ce 
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everywhere,  as  shown  in  Fig.  2.  TJie  notation  Et{Jt,  Me ),  He(Je,  Me)  is  used  to 
denote  the  field  produced  by  Je,  Me  ,  obtained  from  the  usual  potential  integrals. 

In  other  words, 

£«(7«,M«)  =  £e(7e,0)  +  Ee(0,Mj)  (3) 

He(Je,Afe)  =  He(Je,0)  +  He(0>Mt)  (4) 

where 

Ee(7e,0)  =  -ju>Ie(7e)  -  V^(7e)  (5) 

He(7e,0)  =  —  V  X  Ie(7e)  (6) 

Pe 

£«(0,M«)  =  --V  X  Fe(Me)  (7) 

fe 

ffe(0,Me)  =  -jwFe(Mt)  -  V^(A/e)  (8) 

Here  .Ae,  Fe,  ,  and  are  the  usual  magnetic  and  electric  potentials 

Ae  =  Me  JJs  Jt{r)Ge(r  -  r)ds  (9) 

4>e  =  ~  Jfsqe{rl)Ge(r -r)da  (10) 

Fe=ce  JJ  Jte{r)Ge{r  -r'jds'  (11) 

tfe  =  -  //  me^Get?  -  r^)  ds'  (12) 

The  electric  charge  and  magnetic  charge  me  are  related  to  7e  and  Me  by 
the  equations  of  continuity 

qe  =  —  - — V,  •  Je  (13) 

jw 

me  =  -  —  V,  •  Me  (14) 

]U> 

where  V *•  is  the  surface  divergence.  The  Green’s  function  is  that  for  infinite 
media, 

e->*e|r-r'| 

c' = iij?r7r  <15> 


where  ke  =  u^fiece  . 

_ The  boundary  integral  equations  involve  the  tangential  components  of  E  and 

H  over  S.  However,  the  tangential  component  of  H  is  discontinuous  at  a  surface 
current  J ,  and  the  tangential  component  of  E  is  discontinuous  at  a  surface 
current  M .  Hence,  it  is  important  to  evaluate  tangential  components  on  the 
proper  side  of  surface  currents.  For  this,  we  define  a  surface  S+  to  be  just 
outside  5,  and  5—  to  be  just  inside  S  (see  Fig.  2).  The  notation 


4 


Harrington 


denotes  tangential  components  evaluated  on  S-f,  and  the  notation 

hxE;(Je,Mt)  (18) 

n  x  77;(Je,Me)  (19) 

denotes  tangential  components  evaluated  on  5  -  . 

We  also  need  equivalent  surface  currents  (Jj,  A/j)  on  S  radiating  into  media 
having  constitutive  parameters  everywhere,  as  shown  in  Fig.  3.  The  no¬ 

tation  Ei{J i,  Af  j),  H i(J j,  Af ,)  is  used  to  denote  the  field  produced  by  J,,  Afj, 
obtained  from  the  usual  potential  integrals.  In  other  words,  Ei(Ji,Mi)  and 
are  given  by  (3)— ( 15)  with  all  subscripts  e  changed  to  i .  Again  tan¬ 
gential  components  of  E  and  H  are  discontinuous  at  surface  currents  J  and  M 
on  S .  Just  as  in  the  previous  case,  we  define  a  surface  S+  to  be  just  outside  S , 
and  5—  to  be  just  inside  5  (see  Fig.  3).  To  denote  tangential  components  on 
S+  we  use  (16)  and  (17)  with  subscripts  e  replaced  by  i .  To  denote  tangential 
components  on  5—  we  use  (18)  and  (19)  with  subscripts  e  replaced  by  t. 

M<>  £« 

I.U..JT.) 

7  M7.,  At.) 


Figure  2.  The  surface  currents  Je,  Me  radiatingjnto  medium  pe,£e 
everywhere  produce  fields  Ee(Je,Me),  E  e(Je,  Me) . 


n 


Mo 


Figures.  The  surface  currents  J  j,  Af ,  radiating  into  jnedium  m,£i 
everywhere  produce  fields  Ei{  J  j,  Af  i),  H  ,(J ,,  Af  j) . 
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To  solve  the  original  problem  of  Fig.  1,  we  have  to  solve  the  following  mathematical 
problem:  Inregion^e^,  find  the  field  (1)  and  (2),  where  Tmp ,  7Tmp  are  a  known 
field  and  E  ,  77*ca  are  an  unknown  solution  to  MaxweiFs  equations  satisfying 
the  radiation  condition.  In  region  * ,  the  unknown  field  £,  H  is  a  solution  to 
Maxwell’s  equations.  Furthermore,  the  tangential  components  of  E  and  ~H  over 
5+  in  the  external  region  must  equal  the  tangential  components  of  £  and  H 
over  5—  in  the  internal  region. 

By  the  use  of  equivalent  currents  7e,Mt  over  5  radiating  into  media  e  every¬ 
where,  Fig.  2,  we  can  write  a  general  expression  for  the  field  in  region  e  as 


Ee  =  Evnp  +  £e(7e,  Me)  (20) 

¥e=Fmp  +  £e(J<,Afe)  (21) 

By  the  use  of  equivalent  currents  7i,Mj  over  S  radiating  into  media  »  every¬ 
where,  we  can  write  a  general  expression  for  the  field  in  region  i  as 


Ei  =  £j(7j,M,)  (22) 

Hi  =  (23) 

The  formulas  (20)  and  (21)  give  the  field  ~E  ,H  outside  S  (and  on  5+ ),  and  the 
formulas  (22)  and  (23)  give  the  field  E  ,H  inside  5  (and  on  5-).  Continuity 
of  the  tangential  components  of  E  and  H  from  the  external  to  internal  regions 
requires  that 

nx  [£,mp  +  £<+(7e,M<)]  =  nxI,"(Jj,M1)  (24) 

nx  [Fimp  +  7r<+(7e,J?e)]  =  nx  H~ (3 t)  (25) 

Rearranging  (24)  and  (25),  we  have 

n  x  [£e+(7e,Me)  -  £-(7i,A7j)]  =  -n  x  E,mp  (26) 

*  x  [B?(7e,  Me)  -  H~ (7i,  Mi)]  =  -h  X  TTmp  (27) 

These  are  two  equations  to  determine  four  unknowns  (7e,  Me,7,,  Mj) .  We  must 
enforce  two  more  relationships  among  these  four  unknowns  before  a  unique  solu¬ 
tion  can  be  obtained. 

A.  Electric  Current  Formulation 

One  simple  way  to  reduce  the  four  unknowns  to  two  unknowns  is  to  set 

Me  =  Mj  =  0  (28) 

i.e.,  express  the  fields  in  terms  of  surface  electric  currents  only.  Equations  (26) 
and  (27)  then  reduce  to 

nx  [£e+(7e,0)  -  £"(7i,0)]  =  -n  x  £imp 


(29) 
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n  x  |^(7e,0)  -  Hi  (Ji,0)]  =  -n  x  TTmp  (30) 

We  will  show  in  Section  V  that  the  operator  represented  by  (29)  and  (30)  becomes 
singular  at  frequencies  (or  which  5  ,  when  covered  by  a  perfect  electric  conductor 
and  filled  with  the  external  medium,  forms  a  resonator.  Hence,  numerical  solution 
of  (29)  and  (30)  must  fail  in  the  vicinity  of  such  frequencies. 

This  failure  can  be  clearly  seen  from  Fig.  4  of  [4],  where,  in  addition  to  the 
backscattering  cross  section,  the  authors  compute  the  condition  number  of  the 
matrix  of  the  numerical  solution.  The  matrix  becomes  extremely  ill-conditioned 
at  iga  =  2.405 ,  which  is  the  first  internal  resonance  of  an  empty  conducting 
circular  cylinder.  Note  that  the  computation  went  bad  regardless  of  the  losses 
in  the  material  cylinder,  since  it  is  the  external  medium  (in  this  case  free  space) 
which  determines  the  frequencies  of  failure  of  the  equation. 

B.  Magnetic  Current  Formulation 

Another  possible  way  to  reduce  the  four  unknowns  in  (26)  and  (27)  to  two 
unknowns  is  to  set 

7e  =  7j  =  0  (31) 

i.e.,  express  the  fields  in  terms  of  surface  magnetic  currents  only.  This  is  the  dual 
case  to  A  above,  and  the  equations  are  self  dual.  (An  interchange  of  symbols 
according  to  duality  [9,  Sec.  3-2]  produces  equations  of  the  same  mathematical 
form.)  Hence,  the  dual  equations  must  become  singular  at  dual  frequencies,  i.e., 
frequencies  for  which  the  surface  5 ,  covered  by  a  perfect  magnetic  conductor  and 
filled  with  the  external  medium,  forms  a  resonator.  These  resonant  frequencies 
are  the  same  as  those  for  the  electric  current  case,  section  A  above,  and  hence  the 
magnetic  current  formulation  fails  at  precisely  the  same  frequencies  as  does  the 
electric  current  formulation. 

C.  Combined  Current  Formulation 

Another  way  to  reduce  the  number  of  unknowns  in  (26)  and  (27)  to  two  is  to 
set 

7e  =-7i  =  J  (32) 

Me  =  -M,  =  JI  (33) 

In  other  words,  use  a  single  unknown  electric  current  J  to  represent  both  Jt 
and  —  Jj ,  and  a  single  unknown  magnetic  current  M  to  represent  both  A/e  and 
—  Mi.  This  we  call  the  combined  current  formulation.  The  reason  for  choosing 
(32)  and  (33)  is  most  easily  seen  from  the  viewpoint  of  the  field  formation  of 
Section  IV,  where  we  discuss  it  in  more  detail.  It  is  shown  in  Section  1V-A  that 
the  combined  source  formulation  is  equivalent  to  the  combined  field  foruiulation, 
which  has  been  shown  to  be  nonsingular  at  all  frequencies. 

O.  Other  Choices 

Many  other  relationships  among  Je,Me,Ji,Mi  could  be  chosen.  One  choice 
which  would  lead  to  a  formulation  which  has  no  frequencies  of  singularity  is  an 
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extension  of  the  combined  source  formulation  for  conducting  bodies  [lOj.  This 

involves  choosing 

M*  =  aen  x  J t  (34) 

A/j  =  a{h  x  (35) 

where  ae  and  a*  are  constants  to  be  chosen.  For  reasons  discussed  in  [10],  a  good 
choice  is  ae  =  rje  and  Qj  =  rjj ,  where  rje  and  r)i  are  the  intrinsic  impedances 
of  region  e  and  region  t ,  respectively.  Also,  a  proof  similar  to  that  used  in  [10] 
shows  that  the  operator  resulting  from  (34)  and  (35)  has  no  singular  frequencies, 
and  no  mathematical  difficulties  would  be  encountered  in  a  numerical  solution  to 
(26)  and  (27). 

IV.  FIELD  FORMULATIONS 

In  this  formulation  we  use  the  equivalence  principle  [9,  Sec.  3—5)  to  pick  a  given  set 
of  sources  able  to  produce  the  desired  field,  and  then  satisfy  sufficient  boundary 
conditions  on  the  field  to  determine  the  sources.  This  is  the  approach  used  in  [6]. 
It  might  alternatively  be  called  a  direct  application  of  the  equivalence  principle. 

According  to  the  equivalence  principle,  an  electromagnetic  field  can  be  termi¬ 
nated  by  placing  the  required  electric  and  magnetic  surface  currents  on  S .  If  we 
specify  that  the  field  and  sources  outside  5  remain  the  same  as  in  the  original 
problem,  Fig.  1,  then  the  currents 

J  =  nxH  =  hx(Himp +  H'cat)  (36) 

M  =  £xn  =  (Eimp  +  E‘cat)  x  h  (37) 

do  not  change  the  field  outside  5  but  produce  zero  field  internal  to  5 .  Since  the 
field  is  zero  internal  to  5 ,  we  can  change  the  medium  to  any  convenient  value. 
In  particular,  if  we  replace  the  internal  medium  by  the  external  medium,  we  have 
J  and  M  radiating  into  a  medium  with  f*e,Ce  everywhere.  This  gives  us  the 
desirable  situation  that  the  potential  integrals  can  be  used  to  calculate  the  field 
from  J  and  M .  This  procedure  gives  us  the  external  equivalence  of  Fig.  4. 

A  second  application  of  the  equivalence  principle  gives  us  the  internal  equiva¬ 
lence  of  Fig.  5.  We  specify  that  the  field  internal  to  S  remain  the  same  as  in  the 
original  problem,  and  the  field  outside  5  be  zero.  This  requires  the  terminating 
current *  h  x  (-H)  and  (~E)  x  n,  the  minus  sign  resulting  from  the  fact  that 
n  now  points  into  the  region  of  zero  field.  Hence,  the  surface  currents  required 
for  Fig.  5  are  just  the  negative  of  those  for  Fig.  4,  given  by  (3G)  and  (37).  Once 
again  we  can  change  the  medium  in  the  region  of  zero  field  to  any  desired  value. 
In  particular,  we  change  it  to  be  equal  to  the  internal  medium,  so  that  we  have 
— J  and  —M  radiating  into  a  medium  with  /i,,  everywhere.  Hence,  we  can 
again  use  the  potential  integrals  to  calculate  the  field  from  —  J  and  —M. 

We  now  use  the  notation  of  Section  11  to  write  four  boundary  integral  equations 
involving  the  tangential  components  of  E  and  H  on  5—  in  Fig.  4  and  5+  in 
Fig.  5.  Just  inside  S  in  Fig.  4,  we  have  n  x  E  =  0  and  hx  H  =  0,  or 
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n  x  F+(J,M)  =  0  (41) 

n  x  77*  (7,  A/)  =  0  (45) 

Note  that  these  are  four  equations  in  two  unknowns  J,M .  Perhaps  certain 
pairs  of  these  equations  can  be  used  to  compute  J,M .  More  generally,  linear 
combinations  of  them  can  be  used. 

A.  Combined  Field  Formulation 

The  method  that  appears  to  have  been  used  the  most  is  to  reduce  the  set  of 
four  equations  to  two  by  adding  (42)  to  (44)  and  (43)  to  (45).  This  gives  us  the 
pair  of  equations 

n  x  [£7(7,  A/)  +  £*(7,A7)]  =  -n  x  Eimp  (46) 

n  x  [ff7 (7,  M  )  +  77*  (7,  M  )|  =  -»  x  Fmp  (47) 

These  are  very  similar  to  the  equations  obtained  from  the  combined  current  for¬ 
mulation,  Section  III— C,  except  for  the  side  of  5  on  which  the  various  fields  are 
evaluated.  To  be  explicit,  if  we  substitute  (32)  and  (33)  into  (26)  and  (27)  we 
obtain 

n  x  [^(AA7)+^.~(7,M)]  =  -n  x  Tmp  (48) 

n  x  [tf*(7,  A?)  +  £"(7,  M)]  =  -n  x  Tmp  (49) 

We  now  show  that,  in  spite  of  the  different  positions  of  +•  and  -  signs,  (46)  and 
(47)  are  equivalent  to  (48)  and  (49). 

At  any  current  sheet  we  have  [9,  p.  34) 

n  x  [7T(1)  -  7^2)]  =  7  (50) 

[£(l)  -  £<2>]  x  n  =  Af  (51) 

where  E^\  :3  the  field  on  th**  srJeof  S  into  which  h  points,  and  £^2\  77^  ^ 

is  the  field  on  the  other  side  of  S  .  In  terms  of  the  notation  of  this  paper, 

n  x  [fl?(7,  M)  -  777(7,  A/)]  =  7  (52) 

n  x  [Et(J,M)  -  £e'(7, M)]  =  -M  (53) 

Rearranging,  we  have 

nxT77(7,M)  =  nx  777(7,M)  +  7  (54) 

h  x  £e+(7,M)  =  »»  x  £7(7, M)  -  If  (55) 

The  same  equations  are  valid  for  subscripts  e  replaced  by  subscripts  i .  Substi¬ 
tuting  (55)  and  the  same  equation  with  subscripts  e  changed  to  *  into  (48),  we 
obtain  (46).  Substituting  (54)  and  the  same  equation  with  subscripts  e  changed 
to  *  into  (49),  we  obtain  (47).  Hence,  the  combined  field  formulation  of  this 
section  and  the  combined  current  formulation  of  Section  III— C  are  identical. 
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It  is  shown  in  [6]  that  the  combined  field  formulation  gives  a  unique  solution 
at  all  frequencies.  Hence,  the  combined  current  formulation  of  Section  UI-C  also 
gives  a  unique  solution  at  all  frequencies. 

B.  £*fleld  Formulation 

If  we  take  only  the  two  equations  involving  tangential  E ,  (42)  and  (44),  we 
have  the  F- field  formulation 

n  x  E~(7,Af)  =  -n  x  Eimp  (56) 

nxF,+(7,A/)  =  0  (57) 

This  formulation  was  used  in  [5]  for  the  more  general  problem  of  both  dielectric 
and  conducting  cylinders  present.  It  may  seem  strange  at  first  that  this  for¬ 
mulation  makes  no  use  of  the  boundary  conditions  on  tangential  H .  However, 
tangential  Tf  is  determined  by  tangential  E  at  most  frequencies  so  that  (56)  and 
(57)  are  sufficient  to  determine  J  and  M  at  most  frequencies.  However,  we  show 
in  Section  V  that  the  F-field  formulation  will  fail  at  frequencies  for  which  5 , 
when  covered  by  a  perfect  electric  conductor  and  filled  with  the  exterior  medium, 
forms  a  resonant  cavity.  These  are  precisely  the  same  frequencies  at  which  the 
electric  current  formulation,  Section  III— A,  fails. 

C.  H -field  Formulation 

Dual  to  the  F-field  formulation,  we  can  obtain  an  //-field  formulation  by  using 
only  (43)  and  (45)  from  the  set  (42)-(45).  Since  the  //-field  formulation  is  dual  to 
the  F-field  formulation,  it  uniquely  determines  J  and  M  at  most  frequencies.  It 
will,  however,  fail  at  frequencies  for  which  5 ,  when  covered  by  a  perfect  magnetic 
conductor  and  filled  with  the  exterior  medium,  forms  a  resonant  cavity.  Since  the 
field  equations  are  self  dual,  these  frequencies  aie  precisely  the  same  frequencies 
at  which  the  F-field  formulation  fails. 

D.  Other  Choices 

Other  choices  of  two  equations  from  the  set  (42)-(45),  and  linear  combinations 
of  the  equations,  could  be  made.  For  example,  instead  of  the  simple  addition  of 
equations  used  in  Section  IV-A,  we  could  take  the  more  general  linear  combina¬ 
tions 

n  x  [f7(7,M)  +  aEt (7,  Af )]  =  -n  x  F*mp  (58) 

n  x  [77e~(7,Af)  +  /?/7,+  (7,Af)]  =  -n  x  H'mp  (59) 

where  a  and  (3  are  constants  to  be  chosen.  It  is  shown  in  [6]  that  any  choice 
of  a  and  /?  for  which  a/3 *  (*  denotes  conjugate)  is  real  and  po.w:ve  gives  a 
formulation  having  a  unique  solution  at  all  frequencies.  In  particular,  the  choice 

a  =  -liftt  (60) 

P  =  -h/p*  (61) 

gives  the  Muller  formulation  (lj.  This  formulation  has  the  advantage  that  the 
static  electric  field  contribution  to  the  left-hand  side  of  (58)  due  to  the  electric 
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charge  associated  with  J  is  zero.  Similarly,  the  static  magnetic  field  contribution 
to  the  left-hand  side  of  (59)  due  to  the  magnetic  charge  associated  with  M  is  zero. 
According  to  the  last  sentence  on  page  300  of  (lj,  the  singularity  of  the  kernels 
in  (58)  and  (59)  due  to  the  electric  and  magnetic  charges  is  no  more  pronounced 
than  the  reciprocal  of  the  distance  between  the  source  point  and  the  field  point. 
Hence,  the  singularity  that  the  kernels  of  the  integral  equations  (58)  and  (59) 
exhibit  as  the  source  point  passes  through  the  field  point  is  not  as  pronounced 
as  the  singularity  of  the  kernels  of  (46)  and  (47).  Computations  have  shown  [6j 
that  the  use  of  Muller’s  formulation  (this  section)  instead  of  the  combined  field 
formulation  (Section  IV- A)  can  lead  to  more  accurate  solutions  for  low  contrast 
bodies. 

There  are  infinitely  many  other  choices  of  combinations  of  equations  from  the 
set  (42)-(45)  that  could  be  made.  There  may  be  theoretical  and/or  computational 
reasons  for  other  choices,  but  such  reasons  are  not  at  present  known.  We  therefore 
do  not  discuss  any  other  choices  of  combinations  of  equations. 

V.  SINGULARITIES  OF  OPERATORS 

We  stated  in  Section  I1I-A  that  the  electric  current  formulation  failed  at  frequen¬ 
cies  for  which  S ,  when  covered  by  a  perfect  electric  conductor  and  filled  with  the 
external  medium,  formed  a  resonator.  In  Section  IV-B  we  stated  that  the  E-field 
formulation  failed  at  these  same  frequencies.  We  now  prove  these  statements. 

For  these  proofs,  we  use  the  theorem  that,  given  an  equation 

Lf  =  g  (62) 

the  operator  L  is  singular  if  the  corresponding  homogeneous  equation 

Lf  =  0  (63) 

has  a  nontrivial  solution  f0.  For  the  electric  current  formulation,  (29)  and  (30), 
we  must  show  that  there  are  some  frequencies  for  which 

n  x  [Ee+  (7«,0)  -  E~ (7it0)]  =  0  (64) 

n  x  [Ee+(7e,0)  -  H7 (7j,0)]  =  0  (65) 

are  satisfied  by  Je,  J,-  ^0,0.  In  words,  (64)  and  (65)  state  that  we  must  seek 
nonzero  Je  such  that  n  x  Ee  and  n  x  He  are  zero  on  5+ ,  or  nonzero  J j  such 
that  n  x  Ej  and  n  x  Hi  are  zero  on  5— ,  or  both. 

In  the  next  paragraph,  we  assume  that  (64)  and  (65)  are  true.  We  use  (64) 
and  (65)  to  derive  (66)  and  (67).  Although  (66)  and  (67)  are  true,  there  is  no 
assurance  that  they  are  as  restrictive  as  (64)  and  (65).  Therefore,  the  set  of  all 
possible  solutions  to  (64)  and  (65)  will  belong  to  the  set  of  all  possible  solutions 
to  (66)  and  (67). 

Consider  the  composite  situations  where,  in  region  e,  (Ee,(./e,0),  He(Je, 0)) 
exists  in  medium  (nt,ce)  and, in  region  »,  (Ei{Ji,  0),  tfj(Jj,0))  exists  in  medium 
Because  of  (64)  and  (65),  the  tangential  components  of  both  the  electric 
and  magnetic  fields  in  the  composite  situation  are  continuous  across  5 .  Therefore, 
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the  composite  situation  has  no  sources  so  that  its  fields  collapse  to  zero.  As  a 
result,  (64)  and  (65)  imply  that 

n  x  £+(7e,0)  =  0 

_  66) 
n  x  H; (Je.O)  =  0 

n  x  E,- (7j,0)  =  0 

n  x  7/~ (7j,0)  =  0 

Regarding  (66),  refer  to  Fig.  2.  From  our  knowledge  of  cavity  resonators,  we 
know  that  there^are  resonances  within  5,  with  currents  Jt  on  S  such  that 
both  tangential  E  and  tangential  H  are  zero  on  S+ .  Hence,  at  each  resonant 
frequency,  there  will  be  a  non-trivial  Jt  which  satisfies  (66).  A  resonant  frequency 
is  a  frequency  at  which  S ,  when  covered  by  a  perfect  electric  conductor  and  filled 
with  (ji«,e«)  forms  a  resonant  cavity.  Moving  on  to  (67),  refer  to  Fig.  3.  Again 
from  our  knowledge  of  cavity  resonators,  we  know  that  there  can  be  resonances 
within  S .  However,  this  time  we  are  evaluating  n  x  E  and  n  x  ~H  on  5- 
as  denoted  Jby  the  —  superscripts  in  (67).  There  are  no  resonances  for  which 
tangential  E  and  tangential  H  are  both  zero  on  S— .  Hence,  (67)  can  never 
have  a  non-trivial  solution  J  { . 

In  the  previous  paragraph,  it  was  shown  that  the  resonant  frequencies  were 
the  only  frequencies  at  which  (66)  and  (67)  could  have  a  non-trivial  solution 
(7e,7,-).  This  non-trivial  solution  is  called  (Jr,  0)  where  Jr  is  the  resonant 
current.  Obviously,  (J R, 0)  satisfies  (64)  and  (65).  It  turned  out  (66)  and  (67) 
are  as  restrictive  as  (64)  and  (65).  The  solutions  to  (64)  and  (65)  are  the  same 
as  t^*  solutions  to  (66)  and  (67).  Therefore,  the  only  frequencies  for  which  the 
total  operator  of  (64)  and  (65)  is  singular  are  the  resonant  frequencies.  This  is 
the  statement  made  in  Section  III- A. 

The  homogeneous  equations  corresponding  to  the  .E-field  formulation  (56)  and 
(57)  are 

hxEte(J,M)  =  0  (68) 

n  x  E* (J ,M)  =  0  (69) 

In  this  paragraph  and  the  next  two  paragraphs,  we  assume  that  (68)  and  (69)  are 
true  and  investigate  how  (68)  and  (69)  restrict  (J ,M).  Consider  the  composite 
situation  of  Fig.  6  where  (Ej(  —  J ,  —  M ),  Hi(  —  J,—M))  exists  in  (pi,Cj)  internal 
to  S  and  (Ee(J,  M),  He(J,M))  exists  in  (fie,ce)  external  to  _5.  The  field  in 
Fig.  6  is  supported  by  equivalent  electric  and  magnetic  currents  J  and  Ad  on  5 
given  by 

J  =  nx  (H  +  (7,M)-7/~(-7,-M))  (70) 

M  =  -nx  (Ee+(7,A7)-  Er(-7,-A/))  (71) 

In  this  paragraph,  we  find  how  (68)  and  (69)  restrict  J  and  M  of  (70)  and 
(71).  We  have 


7  =  nx  (lft(J,M)- He(J,M)} 
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M  =  -nx  (£e+(J,M)-£e(J,A/))  (73) 

and 

7  =  nx  (74) 

M  =  -n  x  (e,+  (7,A7) -E~(J,M)}  (75) 

Substraction  of  (74)  from  (72)  yields 

J  =  n  x  H?(J ,  M)  +  h  x  H~{J  ,M)  (76) 

where  J  is  given  by  (70).  Subtraction  of  (75)  from  (73)  yields 

M  = -h  x  E+(J ,M)  -  h  x  Ef(J ,~M)  (77) 

where  M.  is  given  by  (71).  Since  there  are  no  external  resonance,  (69)  implies 
that 


n  x  H? (J,  M)  =  0 


Equation  (68)  implies  that 


no  resonance 
resonance 


where  “resonance”  means  that  the  frequency  is  a  resonant  frequency  and  “no 
resonance”  means  that  the  frequency  is  not  a  resonant  frequency.  Substitution  of 
(78)  and  (79)  into  (76)  gives 


j  _  f  0,_  no 
\  <*Jr,  res. 

Substitution  of  (68)  and  (69)  into  (77)  gives 


no  resonance 
resonance 


M  =  0  (81) 

_ In  this  paragraph,  we  express  (J,M)  in  terms  of  the  field  just  inside  5  due  to 

{J , M)  of  (80)  and  (81).  Substitution  of  (78)  and  (69)  into  (74)  and  (75)  gives 

J  =  nxH~(-J,-M)  (82) 

M  =  -n  x  E~ (-J,  -M)  (83) 

Now,  is  the  electromagnetic  field  just  inside  5 

due  to  J  in  Fig.  6,  so  that  (82)  and  (83)  can  be  recast  as 

J  =  hxH~[J,  0)  (84) 

M  =  -n  x  E~(J,0)  (85) 


where  (E^e(j7 , 0 ), /fje(*7 |0j.  **  the  electromagnetic  field  just  inside  S  in  Fig.  6. 
This  field  is  radiated  by  J  of  (80)  on  5  in  the  presence  of  (/i,- ,  e* )  inside  S 
and  (/je,«e)  outside  S.  If  the  frequency  is  a  resonant  frequency,  then  (7,  M)  of 
(84)  and  (85)  is  not  trivial  for  a  ^  0  because  the  non-trivial  J  of  (80)  radiates 
a  non-trivial  field  internal  to  5  in  Fig.  6.  However,  if  the  frequency  is  not  a 
resonant  frequency,  then  J  =  0  so  that  the  composite  field  of  Fig.  6  collapses  to 
zero  because  it  has  no  source.  In  this  J  =  M  =  0 . 
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"E,(7,  J7),2?((7,J7) 


Figure  9.  The  composite  situation  used  to  prove  that  the  £- field  formu¬ 
lation  fails  at  resonant  frequencies. 


VI.  DISCUSSION 

We  have  shown  that  there  are  many  different  boundary  integral  formulations  for 
the  problem  of  electromagnetic  scattering  from  and  transmission  into  a  homoge¬ 
neous  material  body.  One  class  of  formulations,  called  source  formulations,  leads 
to  two  equations  in  four  unknown  currents  7g,Mt,Ji,Mi.  Two  additional  re¬ 
lationships  among  the  currents  must  be  postulated.  When  only  a  single  type  of 
current  is  used,  say  electric  only  or  magnetic  only,  the  formulations  fail  at  certain 
resonant  frequencies.  The  formulation  remains  valid  at  all  frequencies  only  if  both 
types  of  current,  electric  and  magnetic,  are  used. 

A  second  class  of  formulations,  called  field  formulations,  leads  to  four  equations 
in  two  unknown  currents  7,M .  We  must  then  choose  which  two  equations,  or 
which  two  linear  combinations  of  the  four  equations,  to  use.  When  only  two  of 
the  equations  are  used,  say  the  E-field  or  the  E-field  equations,  the  formulations 
fail  at  certain  resonant  frequencies.  These  are  the  same  frequencies  for  which 
the  single-source  formulations  fail.  ’7  hen  we  take  linear  combinations  of  the  four 
equations,  such  that  we  are  in  effect  forcing  all  four  equations  to  be  satisfied,  we 
obtain  formulations  valid  at  all  frequencies. 

An  extension  of  the  source-type  formulation  has  been  made  to  give  a  single- 
equation  formulation  for  the  scattering  problem  [11,12].  This  formulation  is  ba¬ 
sically  an  electric  current  formulation  of  the  type  discussed  in  Section  III-A,  for 
which  an  additional  relationship  has  been  used  to  eliminate  J,- .  The  result  is  a 
single  equation  in  terms  of  7e  to  give  the  scattered  field.  The  solution  does  not 
give  the  transmitted  field  directly.  This  involves  finding  Ji ,  the  source  of  the 
internal  field.  Because  this  single-equation  formulation  is  of  the  electric  current 
type,  it  fails  at  those  frequencies  for  which  5 ,  when  covered  by  a  perfect  electric 
conductor  and  filled  with  the  external  medium,  forms  a  cavity  resonator.  How¬ 
ever,  it  is  possible  to  establish  a  single-equation  formulation  in  terms  of  combined 
sources  to  eliminate  this  problem  [13]. 
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Abstract — The  equivalent  circuit  of  a  via  which  connects  two  senti- 
inflnitely  long  micrustrip  transmission  lines  imbedded  in  a  dielectric 
medium  above  a  ground  plane  is  considered.  The  F-typc  equivalent 
circuit  consists  of  an  excess  capacitance  and  an  excess  inductance.  The 
excess  inductance  of  the  via  is  the  same  as  the  one  computed  in  the  case 
of  free  space  (without  a  substrate).  The  excess  capacitance  of  the  via 
is  computed  quasi-statically  by  the  method  of  moments  and  the  image 
method  from  the  integral  equations.  It  converges  rapidly  as  the  num¬ 
ber  of  image  terms  is  increased.  Parametric  plots  of  the  excess  capac¬ 
itance  of  the  via  are  given  for  reference. 


mlcrosttlp  line  1 


l.  Introduction 

VERY  important  type  of  microstrip  discontinuity  in 
integrated  circuits  is  a  via  connection.  An  often-used 
model  for  analysis  of  this  type  of  discontinuity  is  a  mi¬ 
crostrip  via  connecting  two  semi-infinitely  long  transmis¬ 
sion  lines,  as  shown  in  Fig.  1.  This  problem  has  been 
investigated  in  [  1  )-(3}  in  the  case  of  free  space  (i.e.,  with¬ 
out  a  substrate).  The  via  effect  on  the  transmission  line 
(TEM)  modes  is  approximated  by  an  equivalent  circuit  in 
the  previously  mentioned  literature.  The  components  in 
the  equivalent  circuit  are  computed  by  the  method  of  mo¬ 
ments  quasi-statically.  The  problem  is  more  practical, 
however,  if  the  transmission  lines  and  the  via  are  imbed¬ 
ded  in  a  uniform  dielectric,  as  shown  in  Fig.  2.  Here,  the 
widths  of  the  microstrips  are  assumed  to  be  very  small 
compared  to  the  heights  so  that  the  microstrips  are  re¬ 
placed  by  wires  of  the  equivalent  radii  which  are  one 
fourth  the  widths.  This  equivalence  is  justified  in  (4|.  Let 
/),,  h2  be  the  heights  of  transmission  line  I  (wire  I)  and 
transmission  line  2  (wire  2),  and,  a,,  a:.  a,  be  the  radii 
of  wire  1,  wire  2.  and  the  via  (wire  3),  respectively.  By 
assumption,  /»,,  /m,  /i,  -  h,  »  o i.  «y.  <i\.  The  interface 
of  the  dielectric  of  permittivity  c  is  assumed  to  be  parallel 
to  the  ground  plane.  The  coordinate  system  is  so  chosen 
that  the  x-axis  is  on  the  ground  plane  and  parallel  to  wire 
1  and  wire  2  and  that  the  y-axis  coincides  with  (he  axis  of 
wire  3  (see  Fig.  2). 
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Fig.  I.  A  microstrip  discontinuity:  via  connection. 


Fig  2.  A  transmission  line  model  of  a  microstrip  via. 


The  equivalent  circuit  of  the  via  for  the  quasi-TENl 
modes  is  shown  in  Fig.  3.  where  Zjn  and  Zjn  are  the  char¬ 
acteristic  impedance  of  the  transmission  lines  I  and  2. 
respectively.  Since  the  dielectric  is  nonmagnetic.  L,.  is  the 
same  as  that  in  free  space,  which  was  determined  in  (3|. 
We  wish  to  determine  Cc  in  this  paper. 
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Fig.  3.  The  equivalent  circuit  of  the  via. 

II.  Formulation 

The  excess  capacitance  in  Fig.  3  is  a  quasi-static  quan¬ 
tity  and  is  defined  formally  [I]: 

C,  =  lim  Q  ~  (,) 

/,./>-  ®  v 


where  Q  is  the  total  charge  on  wire  1  of  length  I,  ( wire  I 
truncated  at  x  =  -/, ),  wire  2  of  length  /,  { wire  2  trun¬ 
cated  at  x  =  /,),  and  the  via.  V  is  the  constant  potential 
maintained  on  the  surfaces  of  the  wires  with  respect  to 
ground.  It  is  convenient  to  assume  that  V  is  1  V  through¬ 
out  the  following  text.  <701  and  qn2  are  the  uniform  charge 
densities  of  single  transmission  lines  I  and  2.  Numerical 
computation  of  C,  directly  from  (1)  is  not  recommended 
because  the  numerator  is  the  result  of  subtraction  of  large, 
nearly  equal  numbers  as  /,  and  l,  get  large.  The  more  ef¬ 
ficient  method  is  to  decompose  the  total  charge  density 
into  the  uniform  charge  density  and  the  excess  charge 
density.  The  latter  is  then  determined  by  the  method  of 
moments. 

Let  q  be  the  total  charge  distribution  Then  we  may 
write 


<7m  +  <7,, 

q  =  <7ii2  +  q,. 

V<7„ 


on  wire  I 
on  wire  2 
on  wire  3. 


Here  r/,7 1  is  </ol  existing  on  wire  I  and  qt[ ,  is  qil2  existing 
on  wire  2.  If  q^  is  dclincd  as  </„,  existing  on  the  extension 
of  wire  1  to  x  =  oo  and  if  q„2  defined  as  r/„2  existing  on 
the  extension  of  wire  2  to  x  =  -oo,  we  have 

)  +  'Hqoi )  =  L  on  wire  /.  i  =  1,2 


where  • )  is  the  potential  due  to  the  charge  density 
the  presence  of  the  dielectric  and  the  ground  plane. 

The  total  charge  distribution  q  is  governed  by  the 
lowing  boundary  condition 

i p(q)  =  I ,  on  wires  1,  2.  3. 

In  view  of  (2)  and  (3),  wc  write 

'P(qc)  =  M'/oi)  -  iUtii),  on  wire  I 

<7<- )  =  'j'(qn)  -  <H<7ni).  on  wire  2 

(</,.)  =  1  -  '/'(‘/in)  -  tA ( r/02 ) •  <>n  wire  3. 


Note  that  we  replaced  I  -  ^(<70, )  by  i£(<?oi)  on  wire  1 
and  I  -  i by  ^(q^)  on  wire  2  to  avoid  possible 
significant  round-off  errors. 

Equations  (5)  are  integral  equations  to  determine  the 
excess  charge  distribution.  In  order  to  obtain  a  moment 
method  solution  for  qe,  we  first  truncate  q,  and  (5)  at  x  = 
—  AT|/i,  on  wire  1  and  x  =  K2h2  on  wire  2.  This  approxi¬ 
mation  is  justified  due  to  the  fact  that  qt  and  the  right- 
hand  sides  of  the  first  two  equations  in  (5),  called  the  re¬ 
sidual  potentials  on  wire  1  and  wire  2,  respectively,  decay 
rapidly  along  wire  1  and  wire  2.  Experience  has  shown 
that  K\  and  K2  should  be  between  2  and  3. 

Let 


.V, 

^  All  P/ll* 

on  wire  1 

It  -  l 

iV: 

n  - 1 

on  wire  2 

,Vt 

^  Aijp#ri* 

/i  —  1 

on  wire  3 

where  { /,„,  n  =  1 ,  2,  •  •  •  ,  A/,,  /  =  1 , 2,  3  }  are  expansion 
coefficients  to  be  determined.  p„,  are  pulse  expansion 
functions  defined  as  follows: 

_  f  1  x,„  <  x  <  x„,  +  Ax, 

P"  ^0  elsewhere, 

n  =  I.  2,  •  •  •  ,  N,;  i  =  I.  2 


=  -*,/»■  +  ('!  -  1)  A.r,. 


=  —  a\  +  (  ii  —  1 )  Ax2, 


A.r,  = 


Mi  +  a, 

Nt 


A.r,  = 


K2h2  +  a2 


v„  <  v  <  y„  +  Ay 
elsewhere. 


y„  =  /t,  +  ( n  -  I )  Ay, 


n  =  1.2, 


/i .  -  h-, 


In  essence,  we  divide  the  truncated  wires  1  and  2  and  wire 
3  into  N |,  N2 ,  subsections,  respectively,  and  assume 
uniform  excess  charge  distribution  on  each  subsection.  A 
simple  moment  method  solution  is  to  satisfy  (5)  at  the 
center  of  each  subsection.  This  results  in  the  following 
matrix  equations: 

^\\Q\  +  Lr.Qi  +  Q\  =  P  i 

L2\Q\  +  ^-22  (?2  "h  Li \Q\  =  P2 

L11G1  +  (?2  +  Lyi  Qy  =  Py  (7) 

where  Q,  arc  the  column  vectors  of  the  coefficients  /,„.  L,, 
arc  matrices  whose  mirth  elements  are  the  potentials  due 
to  the  unit  charge  distribution  in  the  nth  subsection  of  wire 
i  at  the  center  of  the  nith  subsection  of  wire  j,  in  the  pres- 
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cncc  of  the  dielectric  and  the  ground  plane.  /*  are  the 
'  column  vectors  whose  with  elements  are  residual  poten¬ 
tials  at  the  center  of  the  wth  subsection  of  wire  j.  Alter 
solving  the  above  matrix  equations,  we  sum  up  the  excess 
charge  and  obtain  the  excess  capacitance  of  the  via.  That 


(*J  Sl.ol 


r7~  »*  0 
'  ground 


C,  =  Z  Ax,  Y  /„,  +  Ay  X  /„,.  (8) 

t  -  l  m  »  \  n  -  t 

In  the  next  section,  we  will  evaluate  and  Pr 
III.  Computation 

To  evaluate  the  moment  matrices  LJS  and  the  excitation 
vectors  Pj,  one  needs  to  find  the  potential  distribution  due 
to  a  uniform  line  charge  in  the  presence  of  the  dielectric 
and  the  ground  plane.  There  are  several  approaches  to  this 
problem.  One  can  place  a  surface  polarization  charge  on 
the  interface  and  solve  the  problem  in  free  space.  This  is 
an  efficient  approach  in  the  two-dimensional  case  (i.e.. 
the  line  is  infinitely  long),  as  (5|  demonstrated.  However, 
in  the  three-dimensional  case  (i.e.,  the  line  is  finite),  one 
needs  a  huge  number  of  expansion  functions  for  the  po¬ 
larization  charge  and  the  method  becomes  too  time  con¬ 
suming.  Another  approach  is  to  use  the  image  method. 
Although  the  image  method  introduces  an  infinite  set  of 
image  line  charges  and  is  difficult  to  generalize  to  the  case 
of  multilayered  dielectrics,  the  solution  is  fairly  straight¬ 
forward  and  is  obtainable  numerically  easily  in  the  case 
of  one  dielectric  interface.  We  take  the  second  approach. 

The  electrostatic  potential  due  to  a  unit  point  charge  in 
the  presence  of  a  dielectric  interface  parallel  to  an  infinite 
ground  plane  (i.e.,  the  Green’s  function  of  our  problem) 


Ui  •Sj.o) 


♦  t0,  *2,0  I 


777  **  0 
'ground 


l-'ip.  4.  Three  cases  ol  line  charge  sources. 

using  an  integral  representation  of  the  free-space  Green’s 
function  |10|.  111).  Then,  the  potential  distribution  due 
to  a  uniform  line  charge  is  obtained  by  integration  over 
the  line  source.  Reference  |6|  gives  the  complete  detailed 
derivation.' 

Case  I:  A  uniform  line  charge  of  finite  length  with  unit 
line  charge  density  lies  on  the  dielectric  interface.  The 
location  of  the  line  is  specified  by  its  two  end  points  ( X,. 
h 0)  and  (X:,  /t,,  0).  as  is  shown  in  Fig.  4(a).  Call  this 
line  charge  X,.  Then  the  potential  due  to  is  given  by 
16j 


i/T  X, )  I  =  — - —  - — —  In 

v  lr,na,r  4 ireQ  I  +  €r 


!{X2  -  xf  +  (  y  -  /»,)*’  +  z”2  +  X,  -  r 


’(X\  —  x)  +  (  y  —  /i|)  +  z*  +  X\  —  x 


2tr  °°,  / 1  -  er\k  —  *)  +  (  v  +  (2A*  +  1  )/i,)  +  z~  +  X2  -  x 

-=-L—  Y  ( - -  )  In  — - - — - - - 

i  +  tr  t-i>  \i  +  «,/  / - ; — ;  "3  ;■ 

V(X,  -  .r)"  +  (y  +  (2k  +  l ) /i, )  +  r  +  X,  -  .< 


*(X.)| 


rinUiclatnc 


4t«0  1  +  er  t  =«  \  I  +  er 


V(X2  -  xf  +  (y  -  (2k  +  l)/»,)2  +  r  +  X:  -  .i 
-  V(X2  -  .r)2  +  (y  +  (2*  +  l)/»,)2  +  r  +  X2  -  j 
(X,  -  xf  +  (  y  +  (2k  +  l)/t,):  +  c3  +  X,  -  ,1 


'(X,  -  .r)2  +  (y  -  (2k  +  !)/»,)'  +  r  +  X,  - 


has  been  obtained  by  the  image  method  in  which  the  im¬ 
age  charges  are  found  directly  |6|,  |8).  by  considering 
flux  lines  rather  than  the  images  themselves  |9|.  and  by 


'Note  that  many  rclcrcm.es  such  as  1 12)  tin  not  give  the  Green's  turn. turn 
in  all  lour  combinations  ol  source  and  lieltl  regions  .Sonic  references  such 
as  | ? |  give  a  complete  l«rm.  which  is  not  the  simplest  lorm. 
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where  r  =  (x,  v.  c )  is  the  field  point  where  the  potential  According  to  (3).  the  potential  d<ie  to  the  charge  density 
is  to  be  evaluated  and  t,  =  e / e,»  is  the  relative  dielectric  r/,H  on  the  line  (  -oo  <  x  £  oo.  y  =  hit  z  =  0)  is  unity, 

constant.  Since  this  potential  is  given  by  the  product  of  <?0i  with  the 

Case  2:  A  uniform  line  charge  of  finite  length  with  unit  right-hand  side  of  (6,  (A-73)]  with  x,  h,  y ,  and  v '  replaced 

line  charge  density  is  placed  inside  the  dielectric.  The  lo¬ 
cation  of  the  line  is  specified  by  its  two  end  points  (X, , 
h2.  0)  and  (X:,  h2.  0),  as  is  shown  in  Fig.  4(b).  Call  this 

line  charge  At.  Then  the  rv^ntial  due  to  is  given  by  _ 

(61 

I  2f  "  /I  - 


•y  A : )  r 


""  4lT€  1 


i  ( Lz±)' 

+  Crl-0  \l  +  «,/ 


In  — 


[.V;  -  xY  +  (  V  -  >U  +  2A/i,r  +  z2  +  X2-x 


X2)l 


'( X2  -  x)~  +  (  V  +  !u  +  2kli,y  +  r  +  X:  -  X 


V(X,  -  xY  +  (  v  +  h2  +  2 khtY  +  ;2  +  X,  -  x 
V<  X,  —  .«)■’  +  (>•  -  ».  +  2Wi,  )■’  +  +  X,  -  x 

.  J.  £  +  •  z 


4tT6  L  >*  0  47T€  A  -  I  \  l  +  €r 


^(X:  -  x)2  +  (y  -  h2  -  2A7t,)‘  +  22  +  X2  -  x 


'U:  -  xY  +  (  y  +  /i,  -  2Jt/i,r  +  z-  +  X,  -  x 


'(X,  -  r)'  +  (  y  +  /i:  -  2A7i|  T  +  r  +  X,  -  x 


'(X,  -  v)‘  +  (  y  -  /»,  -  2A/i,  )*+:*  +  X,  -  x. 


In  (12),  is  the  series  in  (1 1).  by  < 

•  Case  3:  A  uniform  line  charge  of  finite  length  with  unit 
line  charge  density  is  placed  inside  the  dielectric.  The  lo¬ 
cation  of  the  line  is  specified  by  its  two  end  points  (0,  Y,. 

0)  and  (0,  Y2,  0),  as  is  shown  in  Fig.  4(c).  Call  this  line 
charge  X  ,  Then  the  potential  due  to  X  ,  is  given  by  [6j  _ 

*u,)i  =  — £  (— -r 

Irmu.r  4^  |  +  f(.  *  =0  \  1  +  fr 


by  a,,  /»|,  h,.  and  A,,  respectively,  we  obtain 
47re0  I  +  er  [  \  «i  /  t  «  i  \  I  +  t. 


V+i 


/x:  +  ( Y,  -  2kli,  -  yY  +  r  +  Y2  -  2kh,  -  y 


L  V.v2  +  (  +  2A/r,  +  yY  +  c:  +  Y2  +  2kh ,  +  y 

fx~  +  (  K|  +2kh\  +  y)  +  c"  +  K|  +  2Ar/i  (  +  y 
/.v2  +  (  K,  -  2 (t/i,  -  y):  +  c2  +  K,  -  2A7r,  -  y. 


^(x,)|  dl  =r-  i  £  ({ — - 

■  nnJick'ciriw  4ife  t  „  4 7T€  4-1  \  1  +  er 


X2  +  ( K,  +  2A/i  1  -  y)'  +  ;2  +  Y2  +  2kht  -  y 


L  n/.v2  +  ( K:  -  2kh,  +  yY  +  z2  +  Y2  -  2A/i,  +  y 
/x2  +  (  Y,  -  2kli]  +  y)2  +  C:  +  Y,  -  2A7t,  +  y 
/x2  +  ( K,  +  2kli,  -  y):  +  c2  +  Y,  +  2X7/ ,  -  y_ 


In  (14),  Er.i)  is  the  series  in  (13). 
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According  to  (3),  the  potential  due  to  the  charge  density 
•  ^03  on  the  line  (  ~ 00  2S  v  s  00.  y  =  =0)  is  unity. 

Since  this  potential  is  given  by  the  product  of  (/„,  with  the 
right-hand  side  of  16.  (A-75)|  with  .t.  /i,  v.  and  y '  replaced 
by  a2,  A,,  A2,  and  A2,  respectively,  we  obtain 


(16) 


The  two-dimensional  potentials  on  the  left-hand  sides  of 
(15)  and  (16)  are  also  given  by  Weeks  (7). 

Now,  it  is  trivial  to  write  out  the  expressions  for  the 
moment  matrices  and  excitation  vectors.  In  fact,  we  only 
need  to  use  the  formulas  in  the  dielectric  region.  A  small 
special  consideration  must  be  made  when  the  diagonal 
elements  of  L„  are  evaluated.  In  that  case,  the  A  =  0  terms 
in  the  series  in  (10),  (12),  and  (14)  are  replaced  by 
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0  :o  20  30  *0  50 

NumDcr  of  Image  Terms 

f  i^.  3  C  onxcrgcncc  ol  excess  cupacilaiKC  i  <,  ‘  4> 


O'-  -  . . . . . .  V-M 

o  30  oat  3  c  o<  o  as  o  :«  3:7  33#  0  :•  3  10 

K  =  ;*  .  :  1  . *  - M  - 3 

I'ig.  6.  Normalised  excess  capacitance  of  the  wa  t<.  =4.  In  1  2h ,  </,  >  = 
In  (  2ft : /a : )  =  5  ). 

9  =  0 

0[  =  (2k  +  1 )  h{  +  y,'„;  0:  =  (2k  +  l)/i,  -  y'„ 

03  =  2kh,  +  h 2  +  04  =  2kh,  -  l>:  +  y'„, 

05  =  2A/i,  -  It-,  -  y'm;  0„  =  2A/i,  +  /1,  -  y',. 

In  the  above  equations,  the  primed  coordinates  denote  the 
center  of  each  subsection  and  qm  =  <?,„/<  47re„ ).  i  =  1.2. 

Note  that  the  expressions  given  above  are  not  numerically 
most  efficient  but  show  some  patterns  lor  easy  program¬ 
ming.  The  following  section  gives  some  numerical  re¬ 
sults. 


IV.  Ni'Mi  kicai.  Ri.sulis  and  Discussion 

The  computer  program  to  calculate  the  excess  capaci¬ 
tance  of  the  via  was  w  ritten  in  Fortran  and  was  run  in  an 
IBM  AT/PC  machine  with  Microsoft  Fortran.  It  only  took 
10  s  to  run  a  case  (with  29  subsections  and  10  image  terms 
in  the  series).  To  verify  our  code,  we  made  two  basic- 
checks. 

Check  I:  If  the  dielectric  constant  is  I  (i.e..  there  is  no 
dielectric  medium  at  all),  then  the  problem  is  identical  to 
the  one  in  1 3 J .  For  the  following  geometry: 

/r,  =  0.08  m  «,  =  0.002  m 

li:  =  0.04  m  <»i  =  0.001  m 


ih  =  0.002  m 


WAlyO  .1  .</  I  X»  I  SN  <  ACAt  ItANt  I  (It  A  Mil  KIIMKII'  VIA 
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K  =  ^  -  I  i  -  i  - 2  s  - s 

Fig  7  Normalised  excess  capacitance  of  the  via  (<,  =  4.  «,/(/(,  -  /i. ) 

=  005). 


Dielectric  Constant 

-  is  - a  - i  s  - s 


Fig.  8.  Normalised  excess  capacitance  of  the  via  (o, /(/;,  -  /i.)  =  0.05. 
In(2/i,/<i,)  =  in(2/i,/ci, )  =  5) 


we  calculated  the  excess  capacitance: 

C,  =  0.3701  pF 

With  the  same  data,  [3|  gives  C(.  =  0.3776  pF. 

Check  2:  The  excess  capacitance  should  be  convergent 
as  the  number  of  the  image  terms  in  the  series  used  in  the 
calculation  gels  large.  With  the  above  data  and  cr  =  4, 
we  plot  the  excess  capacitance  versus  the  number  of  the 
image  terms  in  Fig.  5.  We  observe  rapid  convergence. 

To  be  consistent  with  |3|,  we  plot  the  normalized  ex¬ 
cess  capacitance  of  the  via  in  various  cases,  as  shown  in 
Figs.  6-8.  Figs.  6  and  7  replace  tigs.  4  and  5  in  [3|  and 
help  design  a  rcllcctionless  via  together  with  the  curves 
of  figs.  6  and  8  in  (3|.  Sec  (3|  for  detailed  discussion. 


In  conclusion,  we  used  a  simple  wire  model  to  compute 
the  excess  capacitance  of  a  via  connecting  two  semi-infi- 
nitcly  long  microstrips.  The  approximation  deteriorates 
when  the  microstrips  are  not  narrow  in  comparison  to  the 
heights  of  the  microstrips  above  ground.  In  this  case,  one 
can  divide  the  microstrips  into  rectangles  or  triangles  and 
can  still  follow  the  same  approach  (the  method  of  mo¬ 
ments  and  the  image  method)  as  that  in  this  paper. 
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A  Stable  Integral  Equation  for  Electromagnetic 
Scattering  from  Homogeneous  Dielectric  Bodies 

JOSEPH  R.  MAUTZ,  SENIOR  MEMF.ER.  IEEE 

Abstract — It  is  shown  that  a  previously  derived  integral  equation  (or 
elect roaugnetk  scattering  from  a  homogeneous  dielectric  body  does  not 
have  a  unique  solution  at  resonant  frequencies  of  the  cavity  formed  by 
making  the  surface  S  of  the  body  perfectly  conducting  and  filling  the 
region  Internal  to  S  with  the  external  medium.  This  integral  equation  was 
formulated  so  that  an  equivalent  electric  current  radiates  in  the  presence 
of  the  homogeneous  external  medium  to  produce  the  scattered  field 
external  to  the  body.  A  combination  of  equivalent  electric  and  magnetic 
currants  to  used  to  formulate  an  integral  equation  whose  solution  is 
always  unique. 

I.  INTRODUCTION 

Consider  electromagnetic  scattering  from  a  homogeneous  dielec¬ 
tric  body  immersed  in  a  homogeneous  medium.  The  surface  of  tb 
body  is  S.  Assuming  that  the  scattered  field  can  be  produce  b. 
equivalent  electric  current  J,  placed  on  S  and  radiating  into  „ 
homogeneous  medium  having  the  material  parameters  of  the  external 
ambient  medium.  Marx  [If  derived  an  integral  equation  for  J,  and 
Glisson  [2]  put  this  equation  into  familiar  notation.  A  resonant 
frequency  is  a  frequency  for  which  S,  when  covered  by  a  perfect 
electric  conductor  and  filled  with  the  external  medium  forms  a 
resonant  cavity.  At  a  resonant  frequency,  this  integral  equation  does 
not  uniquely  determine  the  pure  electric  current  that  supposedly 
produces  the  scattered  field  external  to  S.  Actually,  at  a  resonant 
frequency,  a  pure  electric  current  is  not  sufficient  to  produce  an 
arbitrary  scattered  field;  a  magnetic  current  is  also  needed  [3,  pp. 
351-362]. 

II.  Formulation  of  the  Integral  Equation 

Continuity  of  the  tangential  electric  field  across  5  is  expressed  as 
12,  eq.  (9)) 

Ax(K;  +  E')  =  AxF.,;.  (I) 

Except  for  the  superscripts  plus  and  minus,  our  notation  is  that  of  (2|. 
The  superscript  plus  denotes  field  evaluation  on  the  side  of  5  facing 
the  external  region  and  the  superscript  minus  denotes  field  evaluation 
on  the  side  of  S  facing  the  region  internal  to  the  body.  We  assume  that 
the  electromagnetic  field  (E„  H,)  is  radiated  not  by  J,  alone  but  by 
the  combination  of  the  electric  current  J,  and  the  magnetic  current  an 
x  J,,  both  exactly  on  S: 


E,  =  E,(J,,  afixJ,) 

(2) 

H,=  H,(J,,  afixJ,). 

(3) 
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Here,  a  is  an  arbitrary  constant.  The  source  (J„  aft  x  J,)  is  similar 
to  the  combined  source  in  ]4j.  If  a  =  0,  then  the  integral  equation  in 
the  process  of  being  derived  will  be  the  one  presented  in  |2|  The  field 
<E„  Hrf)  is  radiated  by  the  combination  of  electric  and  magnetic 
currents  J„  and  exactly  on  S: 

=  Ea(J<r,  M,r)  (4) 

Hrf=Ha(Jrf,  Mrf).  (5) 

If,  as  in  (2),  (Jrf,  Mrf)  is  to  radiate  no  field  in  the  external  region,  then 

-ax(h;  +  H1)  (6) 

-(E*  +E')xA.  (7) 

The  original  scattering  problem,  the  simulation  of  the  field  in  the 
exterior  region,  and  the  simulation  of  the  field  in  the  interior  region 
are  as  illustrated  in  [2.  Fig.  I]  with  M,  replaced  by  aA  x  Jr. 

Substituting  (2).  (4),  (6),  and  (7)  into  (1)  and  rearranging  terms, 
we  obtain  an  integral  equation  similar  to  [2,  eq.  (24)|: 

Mrf-AxErf  (J„.  Mrf)=  -Ax|E'  +  E;(AxH\  E'xft)|  (8) 
where 

j„=  -Ax!i;(Jf,  afixJ,)  (9) 

Mrf=AxE;(Jr,  afixJ,).  (10) 

The  homogeneous  equation  associated  with  (8)  is 

Mrf-AxE;(jrfl  Md)  =  0  (II) 

where  3d  and  are  given  by  (9)  and  (10).  If  a  solution  J,  to  (8) 

exists,  then  it  is  unique  if  and  only  if  the  homogeneous  equation  (11) 
with  supporting  equations  (9)  and  (10)  has  only  the  trivial  solution  J, 
=  0. 

If  a  =  0,  then  the  integral  equation  (8)  with  supporting  equations 
(9)  and  (10)  is  the  same  as  (2,  eq.  (24)|  and,  of  course,  the 
homogeneous  equation  (II)  with  supporting  equations  (9)  and  ( l<>)  is 
the  same  as  the  homogeneous  equation  associated  with  (2.  cq.  ( 24 ) | . 
When  the  frequency  is  a  resonant  frequency,  there  is  a  nontrivial 
electric  current  J#  such  that 

«xE;(J»,0)  =  0.  (12) 

fix  H  *  (J*,  0)  =  0.  (13) 

Assuming  that  a  =  0,  it  is  evident  from  (12)  and  (13)  that  J,  =  J*  is 
a  nontrivial  solution  to  the  homogeneous  equation  (II)  with  support¬ 
ing  equations  (9)  and  (10).  Therefore,  even  if  12.  eq.  (24)|  had  a 
solution  at  a  resonant  frequency,  then  this  solution  would  not  be 
unique. 

Can  a  different  choice  of  a  cause  (8)-(  10)  to  have  not  more  than 
one  solution?  Otherwise  stated,  can  a  be  chosen  such  that  the 
homogeneous  equation  (11)  with  supporting  equations  (9)  and  (10) 
has  only  the  trivial  solution  J,  =  0? 

Assuming  that  (9)-(ll)  hold,  consider  the  composite  situation 
where,  in  the  internal  region,  (Erf(Jrf,  Mrf),  Hj(Jrf,  Mr/))  exists  in 
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Gv.  «w)  and,  in  the  external  region,  (£,(J„  ad  x  J,),  H,(J,,ad  x 
J,))  exists  in  <m„  «,)•  The  sources  of  the  above  defined  composite 
field  arc  electric  and  magnetic  currents  j  and  M  on  5  given  by 


J  =  ix(H;<J„aixJt)-H;<J<,.  Pi«rf))  (14) 
I*  -  (E;  ( J„  ak  x  J,) -  E;  ( !*„))  x  ft.  (15) 

Equations  (10)  and  (1 1)  reduce  (IS)  to 

M  =0.  (16) 

Equation  (1 1)  is  recast  as 

nx£;(jrf,  !*„)-<>.  (17) 

Since  there  are  no  external  resonances,  (17)  implies  that 

ftxlljfj,,,  IVI.yHO.  (18) 

From  (18).  we  obtain 

AxHj(3d,  !*„)=  (19) 

Equations  (9)  and  (19)  reduce  (14)  to 

J=0.  (20) 


>  0.  then  the  integral  equation  (8)  with  supporting  equations  (9)  and 
(10)  can  not  have  more  than  one  solution  J,. 
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Surface  Impulse  Response” 
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According  to  (16)  and  (20),  the  composite  field  defined  in  the 
paragraph  containing  (14)  has  no  sources.  Therefore,  this  field 
collapses  to  zero.  In  particular, 

nxE/(J„aAxJ,)  =  0.  (21) 

It  is  shown  in  (4,  eqs.  (4)-(8)|  that,  if  Re(a)  >  0,  then  the  only 
solution  to  (21)  is  J,  ■  0.  Here,  Re  denotes  the  real  pan.  Thus,  if 
Re(a)  >  0,  then  only  J,  =  0  can  satisfy  (9M 11).  Therefore,  if  Re(a) 
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have  appeared  as  a  communication. 
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Abstract  —  A  simple  numerical  procedure  is  implemented  to  find  the  loci  of  the  TE 
and  TM  leaky  wave  poles  for  a  grounded  dielectric  slab  as  the  thickness  is  varied. 
The  information  of  how  these  complex  poles  are  distributed  is  very  important  when 
various  deformed  integration  paths  for  Sommerfeld  integrals  are  considered. 
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1  Introduction 

It  is  well  known  that  the  electromagnetic  fields  due  to  a  dipole  in  a  layered  media 
in  an  open  region  can  be  expressed  in  terms  of  an  improper  integral  [1]  to  account 
for  the  continuous  spectra.  In  order  to  perform  this  integration  accurately  and 
efficiently,  the  integration  path  must  be  deformed  off  the  real  axis  in  some  cases. 
Various  deformed  paths  have  been  considered  by  many  researchers,  e.g.,  Newman 
and  Forrai  [2],  Michalski  and  Zheng  [3],  Fang  and  Chew  [4],  and  Sarkar  [5],  to  name 
a  few.  For  half-space  problems,  there  are  only  two  surface  wave  poles  located  on 
either  Sheets  I  and  IV  or  Sheets  II  and  III  of  the  Riemann  surface  [3,  6],  depending 
on  the  constitutive  parameters  of  the  two  media.  For  a  grounded  dielectric  slab, 
besides  a  finite  number  of  surface  wave  poles,  there  are  an  infinite  number  of  leaky 
wave  poles  [7-9].  If  the  branch  cut  is  properly  chosen  [3,  7,  9],  then  all  the  surface 
wave  poles  are  located  on  the  top  sheet  (proper  sheet  or  Sheet  I),  whereas  all  the 
leaky  wave  poles  are  located  on  the  bottom  sheet  (improper  sheet  or  Sheet  II)  of 
the  Riemann  surface.  The  locations  of  the  leaky  wave  poles  are  immaterial  if  the 
integration  path  stays  on  Sheet  I.  In  order  to  improve  the  computational  efficiency, 
deforming  the  path  to  Sheet  II  is  desirable  in  some  cases,  e.g.,  the  steepest  descent 
method  of  integration  for  far  fields,  where  part  of  the  path  enters  Sheet  II  [9.  10. 
13].  In  this  case,  precise  information  on  how  these  leaky  wave  poles  are  distributed 
is  very  important  unless  the  residues  of  these  leaky  wave  "oles  are  negligible.  Some 
researchers  [10,  11],  when  applying  the  steepest  descent  method  of  integration,  or 
its  approximate  version  the  so  called  saddle  point  method,  only  take  into  account 
surface  wave  poles.  Their  results  are  correct  if  the  dielectric  slab  is  thin  and  if  the 
observation  point  is  not  too  close  to  the  interface,  which  is  shown  in  next  section. 
However,  when  the  thickness  of  the  slab  is  moderate,  or  when  the  observation  point 
is  close  to  the  interface,  some  of  the  leaky  wave  poles,  not  too  far  away  from  the 
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saddle  point,  may  have  been  captured. 

The  locations  of  the  TE  and  TM  leaky  wave  poles  can  be  found  by  solving 
for  the  roots  of  two  simultaneous  transcendental  equations  [7,  9],  which  are  then 
mapped  to  the  complex  plane  of  the  integration  variable.  Although  standard  root 
searching  routines  could  be  applied,  we  gain  no  idea  how  these  roots  are  distributed 
on  the  complex  plane,  since  there  are  infinite  number  of  them.  It  is  the  purpose 
of  this  short  communication  to  show  a  simple  numerical  procedure  of  finding  the 
loci  of  the  leaky  wave  poles  as  the  thickness  is  varied,  which  in  turn  gives  us  the 
desired  information. 


2  Root  Loci 

For  simplicity,  we  assume  that  the  grounded  dielectric  slab  is  non-magnetic  and 
lossless.  The  constitutive  parameters  for  the  half  space  and  the  slab  are  (>u0,  e0) 
and  (/io,e0er),  respectively.  The  thickness  of  the  slab  is  d.  The  locations  of  the 
surface  and  leaky  wave  poles  are  the  roots  of  [9,  10,  12], 

£>te  =  J  l/1  “  £2  +  cot  ( ->"M)  ( 1 ) 

for  the  TE  case,  and 

£>TM  =  ier v/1  -£2  -  \J*r  -£2M)  (2) 

for  the  TM  case.  Here  £  is  the  normalized  parameter  such  that  kp  =  £k0  is  the 
transverse  wave  number  or  integration  variable  for  the  Sommerfeld  integral  in 
terms  of  a  transmission  line  representation  in  the  z  direction  [9,  10],  where  z  is 
perpendicular  to  the  interface.  For  convenience,  the  branch  cut  is  chosen  to  be 
the  line  such  that  Im{ \J\  —  if4}  =  0  [3,  9,  10,  13],  as  shown  in  Fig.  1.  When  this 
branch  cut  is  chosen,  the  roots  of  (1)  and  (2)  on  Sheets  I  and  II  of  complex  £  plane 
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are  the  surface  wave  and  leaky  wave  poles,  respectively.  It  is  well  known  that  the 
surface  wave  poles  are  real  and  confined  in  the  region  (1,  v/c7),  and  the  number  of 
them  is  finite  [10,  12].  Furthermore,  it  can  be  shown  that  the  leaky  wave  poles 
assume  both  real  and  complex  values,  and  the  number  of  them  is  infinite  [7,  9]. 
A  simple  graphical  procedure  [9,  12]  can  be  employed  to  find  the  real  roots  of  (1) 
and  (2),  as  shown  in  Figures  2(a)  and  2(b),  respectively.  In  these  two  figures, 

X  =  k0dyjer  -  <f2  and  B  =  kad\JtT  -  1.  (3) 

The  intersecting  points  on  the  upper  and  the  lower  parts  are  related  to  the  surface 
and  leaky  wave  poles  on  the  Re(£)  axis,  respectively. 

There  is  only  one  pole  associated  with  curve  (1)  in  Fig.  2(a).  It  is  a  surface  wave 
pole  if  B  >  t/2,  a  leaky  wave  pole  in  the  regions  1  <  £  <  y/el  and  y/Tr  <  {  <  oo 
if  1  <  B  <  tt/2  and  B  <  1,  respectively  [7],  [9].  Likewise,  the  pole  associated  with 
curve  (1)  in  Fig.  2(b)  is  always  a  surface  wave  pole.  Clearly,  there  are  two  poles 
associated  with  the  rest  of  the  curves  in  both  figures.  As  an  illustration,  let  curve 
(2)  in  Fig.  2(a)  be  examined.  If  B  >  37r/2,  the  semi  circle  intersects  with  curve  (2) 
at  two  points,  one  for  the  surface  wave  pole  and  the  other  for  the  leaky  wave  pole. 
As  B  is  reduced  and  becomes  smaller  than  Zir/2,  the  surface  wave  pole  falls  down 
to  Sheet  II  from  Sheet  I  through  £  =  I,  which  is  the  branch  point.  It  then  becomes 
a  leaky  wave  pole  and  moves  toward  the  right,  whereas  the  original  leaky  wave 
poles  stays  on  the  same  sheet  and  moves  toward  the  left.  As  B  reaches  a  value  B0, 
such  that  the  semi-circle  is  tangent  to  curve  (2),  these  two  leaky  wave  poles  merge 
into  one.  If  we  reduce  B  further,  then  the  semi-circle  does  not  intersect  with  curve 
(I)  any  more,  and  the  poles  move  symmetrically  off  the  Re(f)  axis  on  Sheet  II.  It 
is  understood  that  the  poles  move  smoothly  as  B  is  varied.  Based  on  the  above 
observation,  the  numerical  procedure  of  finding  the  locus  of  the  leaky  wave  poles 
in  the  £  plane  associated  with  each  curve  in  Figures  2(a)  and  2(b)  is  as  follows. 
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1.  Find  B0  and  the  corresponding  £,  called  £0.  Since  £0  is  real,  this  step  is  straight 
forward. 

2.  Search  for  a  new  root  &  for  the  parameter  Bx  =  B0  —  A B  with  £0  —  j6  as  the 
initial  searching  point,  where  both  A B  and  6  are  small  positive  numbers. 

3.  Search  for  new  roots  £n+i  for  the  parameters  £n+1  =  Bn-ABn,  n  =  1,2,3,-- 
with  as  the  initial  searching  points. 

If  the  A Bn  are  reasonably  small,  a  few  iterations  of  steps  2  and  3  give  satisfactory 
results  using  a  Newton- Raphson  algorithm.  The  same  procedure  can  be  repeated 
to  generate  the  locus  associated  with  each  curve  in  Figures  2(a)  and  2(b).  Since 
when  £  is  a  root  of  (1)  or  (2),  then  —  £  and  ±£*  are  also  roots,  we  only  show 
the  loci  in  the  fourth  quadrant  (these  are  the  leaky  wave  poles  that  matter  for 
a  time  dependence  e-?u'‘).  Figures  3  to  6  show  the  loci  of  the  leaky  wave  poles 
for  tr  =  4  and  9,  respectively.  Careful  scrutiny  of  these  figures  reveals  that  the 
TM  poles  approach  the  real  axis  faster  than  do  the  TE  poles  as  B  increases  from 
zero.  Furthermore,  for  both  the  TE  and  TM  cases,  the  pole  associated  with  the 
lower  numbered  curve  in  Figures  2(a)  and  2(b)  approach  the  real  axis  faster  than 
that  associated  with  the  higher  numbered  curve.  This  explains  why  only  the 
lower  ordered  leaky  poles  are  important  for  a  thin  slab  when  steepest  descent 
path  is  adopted.  A  comparison  of  Fig.  1,  which  shows  steepest  decent  paths  for 
various  observation  angles,  with  Figures  3  to  6  demonstrates  that  many  leaky 
wave  poles  can  be  captured  by  the  deformation,  depending  on  the  observation 
angle  and  physical  parameters  of  the  slab.  The  residues  of  the  leaky  wave  poles 
are  highly  attenuated  in  the  far  field  as  claimed  by  Fang  and  Chow  [4].  In  fact, 
they  completely  ignored  the  residues  of  the  leaky  poles  for  separation  distance 
k0r  >  and  observation  angle  0  =  x/2,  i.e.,  on  the  interface.  However,  it  is  clear 
from  [14]  that  for  some  values  of  the  dielectric  constant  and  height,  the  closeness  of 
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the  leaky  wave  poles  to  the  steepest  descent  path  still  has  to  be  taken  into  account 
even  for  k0r  =  100  and  0  <  0  <  ir/2.  Although  the  physical  parameters  employed 
in  [4]  and  [14]  are  different,  our  purpose  is  to  emphasize  that  care  must  be  taken 
when  ignoring  the  residues  of  the  leaky  waves  poles.  Moreover,  without  knowing 
the  locations  or  the  loci  of  the  leaky  wave  poles,  this  cannot  be  done  satisfactorily. 

3  Concluding  Remarks 

A  simple  numerical  procedure  for  finding  the  loci  of  TE  and  TM  leaky  wave  poles 
as  the  thickness  of  the  slab  is  varied  is  presented.  These  loci  provide  important 
information  when  the  integration  path  of  the  Sommerfeld  integral  for  grounded 
dielectric  slab  problem  is  deformed  into  the  “improper”  sheet  of  the  Riemann 
surface.  The  accuracy  of  the  loci  has  been  checked  extensively  against  contour 
plots  for  expressions  (1)  and  (2)  with  the  B's  as  parameters. 
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Figure  captions 

Fig.  1.  Steepest  descent  paths  in  the  complex  £  plane  for  different  observation  angles, 
where  the  solid  line  is  on  the  top  sheet,  the  dotted  line  is  on  the  bottom  sheet,  and 
the  zigzag  line  is  the  branch  cut. 

Fig.  2.  Graphic  solution  for  the  TE  and  TM  surface  and  leaky  wave  modes. 

Fig.  3.  Loci  of  the  TE  leaky  wave  poles  in  the  complex  £  plane,  tr  =  4. 

Fig.  4.  Loci  of  the  TM  leaky  wave  poles  in  the  complex  £  plane,  er  =  4. 

Fig.  5.  Loci  of  the  TE  leaky  wave  poles  in  the  complex  £  plane,  er  =  9. 

Fig.  6.  Loci  of  the  TM  leaky  wave  poles  in  the  complex  f  plane,  er  =  9. 
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1  Introduction 

The  dispersion  characteristics  of  an  open  microst  been  an  active  re¬ 

search  topic  over  the  past  two  decades.  Various  methods  and  numerical  results 
have  been  reported  [1]— [12].  Although  these  methods  appeared  rigorous,  some  dis¬ 
crepancies  among  their  results  were  noticed  by  Kuester  and  Chang  [6]  in  1979. 
Since  then,  the  accuracies  have  been  improved  a  lot.  However,  the  spectral  do¬ 
main  solutions  are  all  employed  with  entire  domain  expansion  functions  [1],  [2],  [8], 
[12],  whereas  the  space  domain  solutions  are  associated  with  sobdomain  expansion 
functions  [5],  [10],  [11].  In  this  report,  a  spectral  domain  solution  using  nonuniform 
subdomain  basis  function  is  implemented  to  investigate  the  dispersion  character¬ 
istics  of  the  fundamental  mode  of  a  microstrip  line.  Agreement  with  the  existing 
data  in  the  literature  is  demonstrated  in  several  plots,  where  we  have  used  only 
the  longitudinal  current  in  numerical  computation.  This  numerical  scheme  proves 
to  be  an  alternative,  or  rather  a  supplement,  to  various  well  developed  methods  in 
this  area. 


2  Basic  formulation 

The  microstrip  line  under  consideration  is  shown  in  Fig.  1.  The  conductor  sitting 
on  the  air-dielectric  interface  is  assumed  to  be  of  zero  thickness,  and  the  structure 
is  assumed  lossless.  The  ground  plane  and  the  dielectric  slab  are  of  infinite  extent 
in  the  x  and  the  y  direction,  whereas  the  strip  extends  from  — oo  to  oo  in  the 
y  direction.  The  parameters  of  interest  are  the  phase  constant  and  the  current 
denstities  on  the  strip.  For  simplicity,  we  only  consider  the  fundamental  mode, 
which  is  the  lowest  order  mode  propagating  unattenuated  in  the  z  direction.  The 
problem  is  basically  a  two  dimensional  one.  A  phase  factor  e-J/3y,  where  0  is 
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between  the  wave  numbers  of  the  air  and  the  dielectric,  is  common  to  all  the  field 
and  current  distributions,  and  is  the  only  factor  depending  on  the  y  coordinate. 
Hence, 

f(x,y,z)  =  f(x,z)e~}0y  (1) 

follows  through  the  entire  analysis. 

To  account  for  the  continuous  spectrum  of  the  open  structure,  we  express  the 
field  and  the  current  in  terms  of  the  improper  integral  [13,  Ch.  4] 

s(x'z)=hr^  <2> 

which  is  identified  as  an  inverse  Fourier  Transform.  Thus,  the  Fourier  transform 

of  /,  denoted  by  /,  is  obviously  defined  as 

/(£>*)=/  f(x,z)e*xdx  (3) 

J  —OO 


It  is  known  that  the  fields  can  be  decomposed  into  TE  and  TM  components 
with  respect  to  any  direction  [13,  ch.  3j.  The  most  natural  choice  for  our  problem 
is  the  direction  perpendicular  to  the  diele  'ric  interface,  i.e.,  the  z  direction.  The 
z  components  of  the  magnetic  and  electric  vector  potentials,  i.e.,  A(x,y,z)  and 
F(x,y,z),  which  give  rise  to  TM  to  z  and  TE  to  z  fields,  respectively,  can  be 
expressed  by 


A(x,y,z) 


e~j0y  r  B1e-]Ql{z~d)  e~^xd^  d<z 
J—OO 

e~j0y  r  B2CO-*2Z:  e-*rd£,  0  <z<d 
J — oo  sin  ot2d 


(4) 


AOO 

e~10y  /  Dxt-ia^z-d)  e~*xdZ,  d<z 

F{x,y,z)  =  •  rx>  sin  ctoz  e 

e-jpy  /  £>2i - L.  e-**d£,  0<z  < 

k  J—oo  sin  ot.'id 
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where 


£2  +  P2  +  <>i2  =  w2^0e0  ^ 

(2  +  fi2  +  a]  =  u2fi0e0tr 

are  the  dispersion  relations.  In  (4)  and  (5),  B\yB2,D\  and  D2  are  unknown  func¬ 
tions  of  £  that  can  be  determined,  if  desired.  Furthermore,  to  have  field  which 
decays  from  the  dielectric  and  satisfies  the  radiation  condition,  we  require  that 

/?e(ax)  >  0  and  <  0  (7) 

The  electric  Dyadic  Green’s  function  in  spectral  domain,  evaluated  at  the  dielectric 
interface,  has  been  found  by  many  authors,  e.g.,  Farrar  [5],  Fache  [10],  and  Itoh 
[7],  to  name  a  few.  In  spectral  domain,  we  have  [7] 


Zyy(^d)  Zyx{t,d)  Jy(Z,d) 

ZTyti,d)  Zxx((,d)  Jx(t,d) 


Ey(d,d) 

Ex({,d) 


where 


Zyyii,  d )  ^  ^  +  p  Dm  +  Wo  f  +  02  Dg  (9) 

zyx(z,d)  =  zXy{i,d)  =  \—~tz  ~  u^°~d]  (10) 

7  1  q1q2  ,  P2  1 

xxA'  ]  utoe  +  P  Dm  +UJfi°y  +  l3*~D'  (  } 

are  the  spectral  impedance  functions  evaluated  at  the  interface.  In  equations  (9) 

-(H), 

Dm  =  a2  -  jerQi  cot  a2d  (12) 

De  —  ai  —  ja2  cot  a2d  (13) 

The  spatial  counterparts  of  the  spectral  currents  Jy{(,,d)  and  Jx(£,d),  can  be  ex¬ 
pressed  as 

Jy(x)  =  (14) 
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Jx(x)  =  2>n./*n(x)  (15) 

n 

respectively, where  an  and  bn  are  unknown  coefficients.  Using  the  method  of  mo¬ 
ment  [14],  we  enforce  the  boudary  condition 

Kt  —  0>  on  strip  (16) 


by  taking  the  inner  product  of  the  tangential  electric  field  at  the  interface,  obtained 
from  the  inverse  Fourier  transform  of  (8),  with  Jym  and  Jxm  for  all  m.  Doing  so, 
we  obtain  the  homogeneous  matrix  equation 


IMS] 

’  [0]  ‘ 

.  IMS?,] 

(MS) . 

.  ^  . 

.[°1. 

(17) 


The  eigenvalues  could  be  found  by  searching  the  /3’s  such  that  the  determinant  of 
the  matrix  in  (17)  is  zero. 

By  virtue  of  Parsavel’s  theorem 


I  f(x)g‘(x)dx  =  ~  /  fiOs'iOdt. 

J—  oo  Z/T  J — oo 


(18) 


each  element  of  the  matrix  in  (17)  reads 

1  f°°  3. 


MZ  =  -  jT  j;mZpq(t,  d)Jqndt ,  p,q  =  y,x 


(19) 


Notice  that  some  authors  [4],  [8]  have  neglected  the  conjugate  sign  in  (19).  That  is 
because  the  eigenvalues  remain  the  same  if  either  even  or  odd  entire  domain  basis 
functions  are  chosen. 
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3  Numerical  Procedure 


3.1  Integration  path 

Equation  (19)  is  an  improper  integral,  whose  integrand  contains  Qi  and  a2,  both 
being  double-vaiued  functions.  Since  the  integrands  for  p,  q  =  y,x  are  all  even 
functions  of  a2,  only  ax  introduces  the  branch  points  and  the  branch  cuts  to  the  £ 
plane.  For  bound  mode  in  lossless  structure,  the  branch  points  are  always  on  the 
imaginary  axis  of  the  £  plane.  Fig.  2  shows  the  branch  points  and  branch  cuts 
chosen  by  Kowalski  [2],  Fache  [10]  and  Zheng  [11]. 


Since 


Dm  = 


\M2  -  P  -  £2  -  y/fi»  +  £2  -  kl  cot  (y/erk2  -  -  £2  d), 

*2</?2+£2<a02 

-J  V^2  +  P  -  trkl  -  jtrsfp  +  p-h*  COth  (y//P+V-lrk*  d), 

eTkl  <  £2  +  /32 


-jjp  +  e-  kl  -  jV  c  rkl  -P-e  cot  {y/eJl  -  p  -  £2  d), 

De=  _ _ kl  <  02  +  e  <  a2  (2 

-j\/^2  +  Z2  -  kl  ~  jyjp2  -K2  ~  tr kl  coth  (i//32  +  £2  -  e,-fc2  <f), 

tTkl<e  +  P2 

those  roots  of  and  £>e,  corresponding  to  the  surface  I'M  and  TE  waves  of  the 
grounded  dielectric  slab,  respectively,  must  be  in  the  interval  kl  <  d2  +  £2  <  erfc2. 
Letting  [13],  [15] 

X  =  dyJerkl-02  +  £2>O 
A  =  kody/tr  —  1  >  0 
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(22) 

(23) 


the  roots  of 


-tan*  =  (A2-X2),/2  (24) 

and 

-  X  cot  X  =  (A2  -  X2)1/2  (25) 

then  correspond  to  those  of  Dm  and  De,  respectively.  It  is  seen  from  Fig.  3  and 
Fig.  4  that  there  will  be  at  least  one  surface  wave  pole,  which  is  TM0,  on  the  £ 
plane,  whereas  the  cut-off  frequency  for  TE\  mode  is 


f - _ 

U  4dv/eT~"T 


(26) 


To  find  the  desired  phase  constant,  we  must  try  different  0's  in  (17).  Since  the 
phase  constant  of  the  fundamental  bound  mode  is  always  larger  than  that  of  the 
TM0  surface  wave  [16],  [17],  the  surface  wave  poles  are  understood  to  be  on  the 
imaginary  axis  of  the  (  plane.  The  0's  can  be  restricted  in  the  region  between  e* 
and  the  phase  constant  of  the  T M0  mode  in  the  root  searching  process.  Hence,  we 
can  keep  the  path  of  integration  along  the  real  axis.  In  general,  if  the  higher  order 
bound  modes  are  investigated,  and  if  the  frequency  is  greater  than  the  cutoff’s  of 
the  other  surface  wave  modes,  there  are  some  poles  on  the  real  £  axis.  In  this 
case,  the  deformed  path,  as  shown  in  Fig.  5,  is  proposed.  A  major  advantage  of 
deforming  the  contour  off  the  real  axis  is  that  we  do  not  have  to  search  the  surface 
wave  poles.  Slight  inaccuracy  in  the  locations  of  the  poles  on  the  real  axis  may 
result  in  large  errors,  if  real  axis  integration  is  employed.  Of  course,  the  tradeoff  is 
that  the  integrand  becomes  more  complicated  if  we  deform  the  path  to  the  complex 
plane. 
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3.2  Expansion  function 

For  a  single  microstrip  structure,  entire  domain  basis  functions  have  proven  to  be 
efficient.  However,  in  consideration  of  the  extendibility  to  multiple  strip  structure, 
subdomain  basis  functions  may  be  more  suitable.  Furthermore,  to  retrieve  the 
edge  condition  of  the  current,  non-uniform  discretization  of  the  strip  is  employed. 
The  chosen  basis  functions  are 


' 

1,  £n_i  <C  X  <C  Xn 

< 

0,  otherwise 


(27) 


X  «£n— 1 

”  ,  1  <  X  <C  Xn 

Xn  Xn— i 

Jxn(x)  =  "  - — ,  Xn  <  X  <  Xn+i 

xn+l  ~~  xn 

0,  otherwise 


Their  Fourier  transforms  are 


Jm(£)  =  sm(l^/2)e5gxn_1/2 


j  (n  -^1  fsm(£An/2)  jfan_,/a  _  sin  (£An+1/2) 

An  An+1  J 

where  An  =  xn  -  xn_x  and  xn_x/2  =  (x„_i  +  zn)/2.  Thus 


M»  =  Myy  =  -  r  z  F~{m-n) 

Mmn  *  JQ  *VV  c7 


■  *  'V  j  e 

yy  £2 


(28) 

(29) 

(30) 

(31) 


Mv*  -  -Mxy  =  —  f°°  Z  —■ 

Mmn  ~  Mnm  JQ  ^  c3 


F(m,n )  F(m,n  4-  1)]  A  t 

~~KT~  An+i  * 


(32) 
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mn 


—  Mxx  =  —  f°°  Z  — 

lvInm  _  /  “xx  ca 

7T  JO  £4 


F(m,n) 


F(m+l,n)  F(m,n  +  l)^f(m4-l,n+l) 


AmAn  ^m+X^n  AmA 


m  ^n+1 


Am+i  An+1 


di  (33) 


where 


F(m,  n)  =  sin  sin  (  — "  1/2 )  cos  [(xm-i/2  -  *n-i/2)£]  (34) 


3.3  Computations  of  matrix  elements 

It  is  clear  that  the  efficiency  and  the  accuracy  are  based  upon  how  each  matrix 
element,  in  terms  of  the  improper  integral,  is  evaluated.  The  improper  integral  is 
a  function  of  frequency,  dielectric  constant,  width  and  height  of  the  slab,  width 
and  location  of  the  subsegment,  and  0  used  in  the  computation.  The  behavior  of 
the  integrand  is  very  different  from  one  set  of  parameters  to  another.  Hence,  care 
must  be  exercised  in  performing  the  integrals. 

Close  scrutiny  shows  that  all  the  integrands  are  of  the  order  0(1 /£3)  as  <f  — ♦ 
00.  At  £  =  0,  the  integrands  are  of  finite  values.  As  discussed  in  the  previous 
subsection,  if  0  is  restricted  to  be  larger  than  the  phase  constant  of  the  TM0  mode 
of  the  grounded  dielectric  slab,  then  there  is  no  singularity  on  the  real  axis  of  the 
if  plane.  Hence,  the  convergence  of  the  integrals  is  guaranteed,  as  the  physics  of 
the  system  should  be. 

Due  to  nonuniform  discretization  of  the  strip,  Am  ^  A„,  in  general.  The 
oscillatory  terms  F(m,n),  which  is  a  product  of  three  sinusoidal  functions  with 
/  different  periods,  may  behave^  quite  irregular.  The  integration  thus  become  quite 
difficult.  As  a  remedy,  we  decompose  [0, 00]  into  [0,  L]  and  [T,oo],  where  L  is  less 
than  the  first  zero  of  F(m,n).  The  integrand  over  [T,oo]  is  further  decomposed 
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according  to  the  identity 


sin  A  ■  sin  B  •  cos  C  —  -[cos  (A  —  B  —  C) 

+  cos  (A  —  B  4-  C)  —  cos  (A  -f  B  —  C)  —  cos  (A  +  B  +  C )]  (35) 

Each  integrand  now  contains  only  one  sinusoidal  function,  where  its  zeros  are 
periodic.  Therefore,  it  is  easier  to  perform  the  integration.  As  an  example,  for 
each  M™n,  we  have  to  perform  4  integrals  of  the  form 


JL  q 


where  q  is  the  period  determined  by  (35).  We  can  separate  the  above  integral  over 
[Z/,oo]  into  two  parts,  over  [L,U]  and  [[/,  oo],  where 


U  =  max  50 koy/Tr,  J (— )*  +  tTW0  -  & 


Filon’s  algorithm  [18]  is  applied  to  the  integral  over  [£,  U).  For  the  integrals  over 
[U,  oo],  the  asymptotic  expression  of  Zyy(£)  is  substituted  into  the  integrands, 
which  makes  the  integrals  proportional  to 


cos  u 

B\x)  =  /  — r“u 

Jx  U"5 


The  strategy  for  performing  this  integral  depends  on  the  value  of  x.  If  x  <  1,  then 

=  r-3c'w  (39) 


where  C,(x)  is  the  cosin  integral  defined  by 


cos  u  ~ 

Ci(x)  =  /  - du  =  —y-\nx+Y' 

Jx  U  n=l 


(-l)n+1x2n 
2n  •  (2n)! 


where  7  =  0.5772  is  Euler’s  constant.  If  1  <  x  <  307r,  then 

_  ,  ,  ,  rx  1  —  cos  u  . 


,  ,  .  .  [x  I  -  cos  u  , 

/i(x)  =  — 7  —  Inx  +  /  - du 

Jo  u 
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where  the  last  integral  is  performed  using  Filon’s  algothm.  If  x  >  30t,  then 


_  sin x  yk  (  — l)n(2n)?  cos x  ^  (  — l)n(2n  +  1)! 

~  0T  2s  T2n  Ot  2s  _2n+l 

n= 0  x  '£X  n=0 


(42) 


where  N\  and  7V2  are  such  that  (2iV1+2)(2./V1  +  l)  <  x2  and  (2W2-t-3)(2A72  +  2)  <  x2 
respectively.  It  is  known  that  both  series  in  (42)  are  asymptotic  series  which 
diverge  as  N\  and  N2  approach  infinity  [19].  If  x  is  very  large,  3  or  4  may  be 
more  than  appropriate  for  Nx  and  N2.  Furthermore,  using  Filon’s  algorithm  for 
the  integral  over  [L,U],  we  need  to  divide  [L,U ]  into  subintervals,  over  each  the 
Gaussian-Legendre  quadrature  can  be  used.  The  number  of  subintervals  needed 
is  Af  =  [2 (U  —  L)/q ]  -f  1.  If  M  is  greater  than  30,  the  method  of  average  [15]  is 
suggested  to  perform  the  integral  over  [X,  oo].  Usually,  M  is  likely  to  be  large  if 
frequency  is  high,  or  if  the  width  of  the  strip  is  wide.  However,  the  memoth  of 
average  is  not  implemented  in  our  computer  program. 


4  Numerical  results 

A  computer  program  with  only  the  longitudinal  current  has  been  implemented  to 
check  the  validity  and  the  accuracy  of  the  present  numerical  scheme.  The  symmet¬ 
ric  property  of  the  longitudinal  current  has  been  used.  We  restrict  the  numerical 
results  to  some  typical  configurations  that  can  be  found  in  the  literature.  In  Fig¬ 
ures  6  to  11,  the  effetive  dielectric  constants  for  several  structures  are  presented. 
In  these  figures,  the  solid  lines  represent  our  solutions  using  only  5  expansion  func¬ 
tions,  the  long  dashed  line  the  solutions  from  entire  domain  expansion  functions 
of  the  type  of  the  Maxwellian  charge  distribution,  and  the  short  dashed  lines  the 
surface  wave  modes  of  the  grounded  dielectric  slab.  Furthermore,  the  square  and 
triangle  represent  the  results  sampled  from  the  literature.  Attention  could  be  paid 
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to  Fig.  10,  which  structure  has  been  tested  by  many  researchers,  and  theur  dis¬ 
crepancies  were  detected  by  Kuester  and  Chang  [6].  In  Figures  6  and  7,  the  results 
from  subdomain  and  entire  domain  expansion  functions  are  almost  identical,  and 
are  very  close  to  those  obtained  by  Getsinger  [3],  whose  dispersion  model  is  based 
on  approximating  the  microstrip  line  as  suitable  waveguides.  Although  some  de¬ 
viations  between  our  solutions  and  Getsinger’s  are  observed  in  Figures  8  and  9, 
it  does  not  necessarily  imply  that  our  solutions  are  less  accurate.  From  Figures 
10  and  11,  it  is  clear  that  the  subdomain  basis  function  yieldjl  better  results  than 
do^(  tbe  entire  domain  basis  function.  . 

Figures  12  and  13  show  the  normalized  longitudinal  currents  for  d  =  3.048mm, 
w  =  3.175mm  and  er  =  11.7.  In  Fig.  12,  the  square  and  the  triangle  are  the 
currents  at  the  center  at  each  subdomain  for  5  and  10  expansion  functions,  respec¬ 
tively,  whereas  the  solid  line  is  for  20  expansion  functions.  In  Fig.  13,  the  currents 
are  computed  for  three  different  frequencies,  0.001  GHz,  2  GHz  and  12  GHz,  all 
with  10  expansion  functions.  Although  the  frequency  var/ quite  a  lot,  the  current 
profiles  remain  close  to  the  Maxwellian  charge  distribution.  Notice,  we  have  usef 
more  expansion  functions  near  the  edge  than  near  the  center  of  the  strip,  in  order 
to  retrieve  the  edge  condition. 

5  Concluding  remarks 

We  have  shown  that  the  numerical  scheme  based  on  spectral  domain  analysis  with 
non-uniform  subdomain  expansion  functions  could  be  an  alternative,  or  at  least 
a  supplement,  to  the  existing  methods.  Although,  we  only  investigate  single¬ 
strip  structures,  this  approach  is  suitable  to  asymmetirc  multiple-strip  structure? 
as  well.  To  improve  the  accuracy,  the  transverse  current  component  should  not  be 
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neglected.  However,  the  improvement  by  incorporating  the  transverse  current  and 
the  extension  to  multiple  structure  are  warranted. 
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1  Introduction 

Due  to  the  advantage  of  low  cost,  small  size,  light  weight,  and  high  reliability  [1-4], 
the  microwave  printed  circuit  (MPC)  has  become  more  and  more  popular  during 
the  past  three  decades.  The  theories  of  the  MPC  have  been  widely  discussed  in 
the  literature.  An  extensive  listing  of  reference  papers  can  be  found  in  [3].  Most 
models  and  numerical  schemes  that  have  been  proposed  and  published  regarding 
the  MPC  are  valid  at  low  frequencies  [5],  [6].  Many  numerical  techniques,  excluding 
quasi-static  analysis,  are  summarized  in  [7].  A  rigorous  approach  using  an  integral 
equation  formulation,  which  takes  into  account  the  radiation  and  the  surface  wave 
that  are  excited  at  the  discontinuities  of  a  microstrip,  has  been  proposed  by  Katehi 
and  Alexopoulos  [8].  In  this  study,  we  seek  the  same  objectives  but  with  different 
computational  schemes. 

The  basic  structures  under  consideration  are  shown  in  Fig.  1(a)  and  Fig.  1(b) 
for  two  and  three  layered  MPC’s,  respectively.  All  the  media  are  assumed  isotropic, 
lossless,  and  non-magnetic.  The  first  layer  is  free  space,  whereas  the  second  and 
the  third  layers  are  characterized  by  (£i,/i0)  and  (e2 ,Mo),  respectively.  The  planar 
circuits  are  printed  on  the  interfaces.  For  the  three  layered  structure,  a  via  con¬ 
necting  the  circuits  on  the  two  different  interfaces  is  allowed,  although  this  may 
not  be  a  practical  structure. 

The  method  of  moments  [9]  is  employed  to  investigate  various  MPC  problems. 
The  electric  field  integral  equation  is  formulated  in  the  space  domain.  Instead  of 
computing  the  electric  field  directly,  we  use  the  mixed  potential  integral  equation 
(MPIE)  approach  [10],  [11],  where  the  vector  and  the  scalar  potentials  are  com¬ 
puted  in  the  solution  process.  The  mixed  potential  approach  has  been  explored 
very  thoroughly  for  free  space  problems  [9],  and  its  advantage  over  the  original 
EFIE  is  well  known.  Mosig  and  Gardiol  [10]  extended  the  MPIE  approach  to 
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printed  antennas  on  a  grounded  dielectric  slab.  Later.  Michalski  and  Zheng  [11] 
systematically  developed  this  approach  for  problems  of  arbitrarily  shaped  conduct¬ 
ing  bodies  embedded  in  layered  media.  In  general,  there  axe  many  different  choices 
for  the  vector  potentials.  However,  Sommerfeld’s  choice  has  been  found  to  have 
advantages  over  the  others  [11],  and  is  adopted  here. 


2  Green’s  function  and  integration  path 

To  begin  with,  we  only  deal  with  the  two  layered  structure,  as  shown  in  Fig.  1(a). 
Due’  to  an  x  directed  electric  dipole  on  the  interface,  various  components  of  the 
magnetic  vector  potential  and  the  scalar  potential  are  as  follows  [10],  [11],  [16]. 


r 


exp  {-jkzoz) 
sin  kzi(z  +  d )/  sin  kzid 

exp  {-jkzoz) 
cos  kzi(z  +  d)/  cos  kzld 


dkp, 


dkp, 


(1) 


(2) 


G*  oc  /  J0(kpp ) 

Jo 


kP\jkzo  -  kgi  tan  kzld\  j  exp  ( -j kzoz) 

DeDm  1  sin  kzi(z  +  d)/ sin  kzld 


dkp ,  (3) 


where  the  upper  and  the  lower  terms  in  the  curly  brackets  are  for  z  >  0  and 
— d  <  z  <  0,  respectively.  Furthermore,  the  roots  of 


Dz  — -  j kZQ  -f-  kz\  cot  kz\d 


(4) 


and 

Dm  =  je.-kz0  -  kzl  tan  kzld,  (5) 

respectively,  correspond  to  the  TE  and  TM  surface  or  leaky  wave  poles  in  the 
complex  kp  plane.  Notice  that  ±k0  are  the  only  branch  points.  In  general,  the 
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branch  cut  is  artificial  and  can  be  chosen  arbitrarily.  However,  the  most  natural  one 
is  shown  in  Figures  2  and  3,  in  which  the  numbers  of  the  surface  wave  poles,  located 
on  the  upper  half  sheet  of  the  Riemann  surface,  and  the  leaky  wave  poles,  located 
on  the  bottom  sheet  of  the  Riemann  surface,  are  finite  and  infinite,  respectively 
[12],  [13].  If  the  integration  path  is  deformed  off  the  real  axis,  the  residues  of 
the  poles  captured  must  be  taken  into  account.  Some  researchers  [10],  [14]  have 
implemented  the  real  axis  integration  algorithm,  which  is  not  easy  to  perform  in 
two  extremely  cases,  i.e.,  if  B  =  k0dy/er  —  1  is  very  small  or  very  large.  If  the 
parameter  B  is  very  small,  the  first  TM  surface  wave  pole  is  very  close  to  the 
branch  point  kQ,  where  the  integrand  exhibits  a  derivative  discontinuity  [10],  [16]. 
If  B  is  very  large,  then  some  leaky  wave  poles  are  close  to  the  real  axis  [15],  malting 
the  integrand  oscillate  irregularly.  We  choose  the  integration  path  shown  in  Fig.  2 
to  avoid  the  complication  caused  by  the  poles.  However,  if  kDp  gets  larger,  this 
integration  scheme  gets  less  efficient  due  to  the  oscillation  of  the  Bessel  function 
along  the  path.  For  kQp  >  5,  we  make  the  following  substitution  for  the  integrands 
of  (1)  to  (3), 

UKp)  =  \  -  WH-Kp)}  (6) 

The  integration  path  is  then  deformed  to  be  along  the  branch  cut  as  shown  in 
Fig.  3.  Choosing  this  path,  Mosig  and  Gardiol  [16]  has  shown  that  the  integral 
can  be  decomposed  into  three  parts,  i.e.,  the  static  term  which  is  purely  real,  the 
space  wave  term,  and  the  surface  wave  term.  If  k0p  is  very  large,  the  efficiency 
is  again  degraded  by  the  rapid  oscillation  of  the  Bessel  function  in  the  interval 
kp  6  [0,  lb0].  In  this  case,  the  steepest  descent  path  may  be  a  better  candidate. 
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3  Numerical  results  for  printed  dipoles 

Figures  4  to  9  show  the  current  distribution,  and  the  input  and  mutual  impedances 
of  some  printed  wire  and  strip  dipoles.  The  results  agree  with  those  available  in 
the  literature  [13],  [17],  [18].  Attention  should  be  paid  to  Fig.  10,  where  we  use  49 
unknowns  in  the  solution  process  for  a  strip  dipole  of  width  0.005Ao,  and  Barkeshli 
uses  only  three  unknowns  for  a  dipole  of  width  0.01Ao.  Inaccuracy  in  Barkeshli’s 
solution  is  suspected,  since  in  Fig.  10(a),  the  results  can  hardly  be  differentiated 
from  those  of  a  wire  dipole  of  radius  0.00125Ao. 

4  Proposed  study  for  microwave 
printed  circuits 

In  the  MPC,  a  microstrip  line  always  serves  as  a  transmission  line.  The  number 
of  ideal  transmission  lines  needed  to  represent  the  microstrip  must  be  the  same 
as  that  of  the  bound  modes  that  can  be  excited  [3].  For  simplicity,  we  restricted 
ourselves  to  the  frequency  range  in  which  the  fundamental  mode  is  the  only  bound 
mode.  Furthermore,  we  assume  that  w  <<  A0,  where  w  and  A0  are  the  width  of 
the  microstrip  and  the  free  space  wavelength,  respectively.  Hence,  the  transverse 
current  can  be  ignored  with  little  error  introduced.  The  axial  current  has  singu¬ 
larities  at  the  two  edges  of  the  microstrip.  The  transverse  variations  of  the  axial 
current  at  several  frequencies  are  shown  in  Fig.  11,  where  the  Maxwellian  function 
is  plotted  for  comparison.  In  general,  if  the  frequency  gets  higher,  the  current 
distribution  deviates  from  the  Maxwellian  function  more.  Figures  12  and  13  give 
curves  of  (A<,/Aj)2  versus  frequency,  which  are  compared  with  those  obtained  by 
Getsinger  [19],  Farr  air  [20],  and  Fache  and  De  Zutter  [21].  It  is  obvious  that  if 
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w/d  <  1  and  if  the  frequency  is  not  too  high,  the  Maxwellian  function  for  the 
current  gives  a  reasonably  accurate  result  for  the  propagation  constant.  To  further 
improve  the  accuracy,  we  assume  the  transverse  variation  of  the  current  to  be  of 
the  form 

J'(y)  CC  1  [T.( 2y/w)  +  W, (2y/w))  (7) 

V 1  -  Psi/®)2 

where  Tn  is  the  Chebyshev  polynomials  of  the  first  kind  of  order  n.  After  the 
propagation  constant  is  found,  the  unknown  coefficient  b  can  be  easily  solved.  The 
characteristic  impedance  can  be  computed  according  to  [22] 


where  P  is  the  real  power  propagating  in  the  axial  direction.  Fig.  14  shows  the 
normalized  propagation  constant  and  the  characteristic  impedance  computed  from 
one  (b  =  0)  and  two  (b  0)  expansion  functions.  The  results  computed  from  the 
commercial  software  SUPERCOMPACT  and  from  quasi-static  analysis  are  also 
plotted  for  comparison.  The  power  ratio  in  the  dielectric  and  the  b  in  (7)  are 
shown  in  Fig.  15.  The  modeling  of  the  discontinuity  is  illustrated  in  Figs.  16 
and  17.  By  observing  the  standing  wave  pattern  exhibited  in  Fig.  17(a),  several 
parameters  can  be  computed.  Notice  that  we  need  to  compute  two  sets  of  currents 
with  two  different  Li s  in  Fig.  17(a)  in  order  to  determine  the  scattering  matrix. 
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(a)  Two  layered  configuration 


(W 
d,  <VV 
d2  (e2’  V 


(b)  Three  layered  configuration 


Fig.  1.  Side  view  of  microwave 
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circuit. 


Fig.  2..  Deformed  integration  path  in  the  complex  kp  plane. 


Fig.  3.  Integration  path  along  the  branch  cut  in  the  complex  k.  plane. 


Current  distribution  of  a  center-fed  wire  dipole,  L  =  0.5Ao,  er 
0.1016Ao,a  =  0.0001  A0,  (a)  our  result,  (b)after  Rana  (1981). 


Current  distribution  of  a  center-fed  wire  dipole,  L  —  0.5Ao,£ 
0.15Ao,a  =  0.00005Ao,  (a)  our  result,  (b)after  Rana  (1981). 


Mutual  impedance  between  two  collinear  wire  dipoles, 
0.127Ao,a  =  0.000lAo,  L  =  0.4Ao,  (a)  our  result, 
(b)after  Alexopoulos  (1981). 


ME 


Mutual  impedance  between  two  broadside  wire  dipoles,  eT 
0.1016Ao,a  =  0.0001  A0,  l  =  0.333A,,  (a)  our  result, 
(b)after  Alexopoulos  (1981). 


Fig.  10 L  Input  impedance  of  the  center-fed  strip  and  wire  dipoles,  er  =  3.25,  d 
0.0796Ao,  w  =  0.005A„,  a  ~  0.00125Ao,  (a)  our  result  with  strip  ( — 
and  wire  ( - ),  (b)after  Barkeshli  (1988). 
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damental  mode  for  a  microstrip. 


Dispersion  characteristics  of  the  fundamental  mode,  where  A0  and  Af 
are  the  free  space  and  guide  wavelengths,  respectively. 


Dispersion  characteristics  of  the  fundamental  mode,  where  XQ  and  A, 
are  the  free  space  and  guide  wavelengths,  respectively. 
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SUPERCOMPACT 
Two  expansion  functions 
One  expansion  function 
Quasi-static  analysis 


SUPERCOMPACT 
Two  expansion  functions 
One  expansion  function 
Quasi-static  analysis 


Fig.  f4  (a)  Normalized  propagation  constant  and  (b)  characteristic  impedance 
for  a  microstrip  line  with  e,  =  9.6,  w  =  0.6mm,  d  =  0.6mm. 
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(a)  subdomain  basis  (unction  (or  smaN 
double  step  discontinuity 
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(b)  subdomain  basis  function  for  large 
double  step  discontinuity 


<c)  equivalent  network  represented  by  its 
scattering  matrix 

Fig.  l€  Modeling  of  microstrip  double  step  discontinuity. 
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